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Summary. We review the role of strong electronic correlations in quasi-two- 
dimensional organic charge transfer salts such as (BEDT-TTF)2X, (BETS)2^ and 
/3'-[Pd(dmit)2]2-Z- We begin by defining minimal models for these materials. It is 
necessary to identify two classes of material: the first class is strongly dimerised 
and is described by a half-filled Hubbard model; the second class is not strongly 
dimerised and is described by a quarter filled extended Hubbard model. We argue 
that these models capture the essential physics of these materials. We explore the 
phase diagram of the half-filled quasi-two-dimensional organic charge transfer salts, 
focusing on the metallic and superconducting phases. We review work showing that 
the metallic phase, which has both Fermi liquid and 'bad metal' regimes, is described 
both quantitatively and qualitatively by dynamical mean field theory (DMFT). The 
phenomenology of the superconducting state is still a matter of contention. We criti- 
cally review the experimental situation, focusing on the key experimental results that 
may distinguish between rival theories of superconductivity, particularly probes of 
the pairing symmetry and measurements of the superfluid stiffness. We then discuss 
some strongly correlated theories of superconductivity, in particular, the resonat- 
ing valence bond (RVB) theory of superconductivity. We conclude by discussing 
some of the major challenges currently facing the field. These include: parameteris- 
ing minimal models; the evidence for a pseudogap from nuclear magnetic resonance 
(NMR) experiments; superconductors with low critical temperatures and extremely 
small superfluid stiffnesses; the possible spin-liquid states in k-(ET)2Cu2(CN)3 and 
/3'-[Pd(dmit)2]2-Z; and the need for high quality large single crystals. 



1.1 Motivation and scope 

The Drudc, Sommerfeld and Bloch models explain many of the properties of 
elemental metals and alloys in terms of a theory of non-interacting electrons. 
This is somewhat surprising as for typical metals the radius of a sphere whose 
volume is equal to the volume per conduction electron is of order 1 A [1] . Thus, 
one naively expects that the Coulomb interaction between electrons will be 
large. Therefore, one of the most significant achievements in condensed matter 
theory is Landau's theory of Fermi liquids [2] which explains, through the idea 
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of adiabatic continuity, why the Coulomb interation simply renormalises the 
weakly interacting electron gas. 

However, in many "strongly correlated" systems interactions qualitatively 
alter the behaviour of the material. For example, in the Kondo effect the inter- 
actions between the conduction electrons and magnetic impurities, which leads 
to a logarithmic increase in the resistivity at low temperatures, cannot be de- 
scribed by perturbation theory from the Fermi liquid state [3] . Indeed even the 
BCS state [4] is not adiabatically connected to the Fermi liquid ground state. 
The Kondo effect is now well understood, but many other strongly correlated 
phases are not, e.g., high temperature superconductivity and the pseudogap 
in the cuprates [5] , colossal magnetoresistance in the manganites [6] and emer- 
gence of unconventional superconductivity in heavy fermion materials [7]. 

The non-perturbative nature of strongly correlated materials presents sev- 
eral major challenges to theory. For example some big questions are: What 
is the appropriate description of the metallic state when the system is near 
to an instability to a Mott insulating state? What are the connections be- 
tween magnetic ordering and unconventional superconductivity? Under what 
conditions can an antifcrromagnct have a spin- liquid ground state and/or 
spinon excitations? When superconductivity is found near a Mott transition, 
what is the relationship between the symmetry of the superconducting state 
and the ground state of the parent Mott insulator? What is the appropri- 
ate microscopic description of a metallic state with a pseudogap? What are 
the mechanisms responsible for unconventional superconductivity and what 
are the appropriate microscopic descriptions of such states? Many of these 
questions are deeply intertwined. 

Organic charge transfer salts are excellent model systems in which to 
study many of the questions about the strongly correlated phenomena de- 
scribed above. Band structure suggests that these materials should be metals 
at all experimentally relevant pressures and temperatures. However, as we 
will discuss below, the observed phases include, Mott insulators, Neel anti- 
ferromagnets, spin-liquids, (unconventional) superconductors, Fermi liquids, 
a pseudogap and a 'bad metal'. Further, organic chemistry allows the or- 
ganic charge transfer salts to be subtly tuned in ways that have never been 
achieved in inorganic materials. A rather beautiful example of this is that in k- 
(ET)2Cu[N(CN)2]Br the Mott transition can be driven by replacing the eight 
hydrogen atoms in the ET molecule with deuterium [8] . Indeed the chemistry 
can be controlled to such an extent that the number of H/D atoms on each 
ET molecule can be varied uniformly throughout the entire sample allowing 
one to move gradually across the Mott transition and observe the coexistence 
of the insulating and superconducting phases expected because the transition 
is first order [9], 

We will argue below that, despite the chemical complexity of the organic 
charge transfer salts, the physics boils down to that of the Hubbard model 
on various lattices and at either half or one quarter filling. We will show 
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that the strongly correlated physics of the various polymorphs and chemical 
constituents can be described within this framework. 
The Hamiltonian of the one-band Hubbard model is 

^Hubbard = - 53 + M%)cL<V + U 53 fl^fl^, (1.1) 

where, tij is the hopping integral from the site i to the site j, \i is the chemi- 
cal potential, U is the Coulomb repulsion caused by putting two electrons on 
the same site, c^J annihilates (creates) an electron on site i with spin a and 
the number operator hi a — c^Cja-. The Hubbard model is, perhaps, the sim- 
plest model that can describe strongly correlated physics and is therefore an 
important starting point for a complete and general description of strong cor- 
relations. Further, the Hubbard model is also believed to describe the essential 
physics of, for example, the cuprates [10] and the cobaltates [11]. 

At some points later in this review we will also wish to discuss extended 
Hubbard models. In particular we will discuss the Coulomb repulsion when 
two electrons are placed on neighbouring sites, we will refer to these as 'V 
terms and the additional term in the Hamiltonian will be 

Hv = V^2(hif +nu)(h n +n jl ), (1.2) 

where the angled brackets indicte the sum is over some specified set of neigh- 
bouring sites, e.g., nearest neighbours only. Another important term we may 
wish to include is the Hcisenberg exchange or 'J' term 

Hj = J]T S. t -S„ (1.3) 

m 

where the spin operator is = J2 a p ^la^apCi^ and er = (a x , a y , o~ z ) is the 
vector of Pauli matrices. Lastly we may wish to include phonons and the 
clcctron-phonon interaction in the model. The simplest approach is to include 
the 'Holstcin' terms which treat dissipationlcss phononic modes via the terms 

^Holstoin = 53 9v («L + O-iv) ^Ha + 53 ^"L^> (I- 4 ) 
iav iv 

where the a^J operators annihilate (create) a phonon in the v th mode on the 
i th lattice site, w„ is the characteristic frequency of v mode and g v is the 
electron-phonon coupling between an electron on site i and a phonon in the 
v th mode on site i. 

This review is not intended to be comprehensive. Rather, we present our 
views of the physics relevant to understand these materials. This is not to say 
that we are the originators of all of these ideas, but rather to stress that, as 
with any healthy research field, there arc a number of prominent researchers 
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who might disagree with the views expressed here. Space constraints will not 
permit us to discuss rival theories and interpretations of data at length, we 
therefore refer interested readers to some alternative points of view [12-14] 
and some more comprehensive reviews [15]. 

1.2 Minimal models 

In this review we will mostly consider materials with chemical formulae of the 
form D%X where D is an organic donor molecule, for example bis(ethylene- 
dithio)tetrathiafulvalene [Ci S 8 H 8 ] (often abbreviated as BEDT-TTF or ET) 
or bis(ethylenedithio)tetraselenafulvalene [CioSegHg] (BETS), and X is an 
anion, e.g., I3 or Cu[N(CN)2]Br. The anion accepts one electron from a pair 
of donor molecules which leads, at the level of band theory, to an insulat- 
ing anionic layer and a metallic donor layer. The chemical complexity of the 
organic charge transfer salts mean that few first principles calculations have 
been reported [16-20]. However, these suggest that a tight-binding or Huckcl 
description of the band structure is a reasonable approximation. In this tight 
binding description a molecule serves as a 'site' [21]. 

1.2.1 Quarter filled charge transfer salts - the (3" and phases 

It is clear that in the model proposed by Kino and Fukyama [21] (described 
above) one expects that there should be, on average, half a hole per site. 
This is indeed the case for the (3" and 9 polymorphs. These show a subtle 
competition between metallic, charge ordered insulating and superconducting 
phases [22,23]. It has been argued that this arises due to the competition of 
spin and charge fluctuations. These effects have been studied in the extended 
Hubbard (t-J'-V-U) 1 model on a square lattice [22,23]. The importance, for 
the superconducting state, of the fact that the actual lattice has a rather lower 
symmetry than the square lattice has also been stressed recently [24] . 

1.2.2 Half-filled charge transfer salts - the (3, /3', k and A phases 

In the /3, /?', k and A polymorphs there is a single intcrmolccular hopping inte- 
gral that is significantly larger than the others. This led Kino and Fukuyama 
[21] to propose that these two molecules behave as a single site as the hopping 
between them is rapid enough that it can be integrated out of the effective 
low energy theory. The two molecules are often referred to as a 'dimer'. How- 
ever, it is important to note that this does not imply that the molecules are 
covalently bonded in the conventional use of the word; i.e., there are no C-C, 

x The J' indicates that the exchange term is to next nearest neighbours. The J 
term can be neglected because in (or near) the charge ordered phase the are no (few) 
occupied sites with occupied nearest neighbours. 
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Fig. 1.1. A sketch illustrating why the band structure of the half- filled layered or- 
ganic charge transfer salts is approximately that of the anisotropic triangular lattice. 
Black circles indicate the dimers, the hopping integral between the two molecules 
in a dimer is the largest energy scale in the problem and can be integrated out of 
the effective low energy Hamiltonian leaving the dimers as effective sites in a Hub- 
bard model description of the electronic degrees of freedom. The largest hopping 
integrals between the dimers, t and t' , are indicated by the maroon and blue lines. 
These form an anisotropic triangular lattice. In this review we take the x and y axes 
of the anisotropic triangular lattice to lie along the directions of the two t hopping 
integrals as this is conventional in the field. This figure is based on the crystal- 
lographic data of Rahal el al. for k>(ET)2Cu(NCS)2 (a typical half-filled organic 
charge transfer salt) [26]. 

S-S or S-C bonds between the two molecules. Rather, the 'bond' is between 
the molecular orbitals themselves, this 'bond' results because of the signifi- 
cant overlap between the highest occupied molecular orbitals (HOMOs) on 
the two molecules, and the fact that the HOMOs are partially occupied. Thus 
the 'bond' between the two molecules forming the dimer is highly analogous 
to a covalcnt bond between atoms. The dimer structure is sketched in figure 
1.1, also shown are the interdimer hopping integrals. Kino and Fukuyama 
gave a parameterisation of t and t' , the interdimer hopping integrals, in terms 
of the intermolecular hopping integrals. It was later realised that this needed 
to be corrected to allow the strong Coulomb repulsion between two holes on 
the same molecule [25]. The interdimer hopping integrals form an anisotropic 
triangular lattice. Thus geometrical frustration may be expected to play a 
significant role if t' ~ t. 

The phase diagram of the k phase materials (the best studied of the half- 
filled organic charge transfer salts) is sketched in Fig 1.2. It is clear that 
this phase diagram cannot result from the non-interacting model described 
above. However, there is a strong Coulombic repulsion between two holes on 
dimer [27]. Therefore, we believe that the Hubbard model may be the simplest 
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model which contains the physics required to describe these systems [21,25]. 
Kanoda [28] proposed that the role of pressure is to decrease U/W where W is 
the bandwidth. As different anions lead to different behaviours (and different 
unit cell volumes) in much the same way as pressure, this can be thought of 
as a effective 'chemical pressure'. This idea has proven an extremely powerful 
framework in which to assimilate the great deal of experimental data now 
available for these systems. 

An often used estimate of U comes from modelling the dimer in a two site 
Hubbard model with an intermolecular hopping integral, t m , and a Hubbard 
repulsion U m for having two holes on the same molecule. In the limit U m 3> 
4t m this model yields U ~ 2t m [25]. However, this is rather misleading because 
quantum chemistry calculations suggest that the Coulomb repulsion between 
holes on neighbouring molecules in the same dimer V m is of the order of U m 
[27]. Thus this term should be included and therefore the extended (t m -U m - 
V m ) Hubbard model should be used to estimate U [29]. It is straightforward 
to show [29] that for this model yields 

U = - (u m + V m - [{U m - V m } 2 + 164] V2 + 4*m) (1-5) 

and in the appropriate limit, (U m + V m ) 3> (U m — V m ),t m , we obtain 

U « i (U m + V m ) > 2t m . (1.6) 

Taking the values of U m , V m and t m reported by Fortunelli et al. [16] we see 
that equation (1.5) gives U = 3.7 eV, this is approximated reasonably well by 
equation (1.6) which gives U = 4.8 eV whereas 2t m = 0.56 eV which is the 
wrong order of magnitude. Thus we see that the plain two site Hubbard (t m - 
U m ) model significantly underestimates U. A more detailed discussion of the 
calculation of the parameters of the Hubbard model will be given in Sect. 1.8.1. 
In particular that section will discuss the fact that the U in vacuo (equation 
1.5) is significantly larger than that in the crystal due to the polarisability of 
the lattice. 

The majority of the remainder of this review with be devoted to these 
half-filled systems. But first we briefly mention three other classes of organic 
charge transfer salts. 

1.2.3 Half-filled charge transfer salts with magnetic anions 

An interesting class of organic charge transfer salts have been prepared 
with magnetic anions. Prominent examples include K-(BETS)2FeCi4 and A- 
(BETS^FeCLt. Perhaps the most novel phenomena observed in these salts is 
field induced superconductivity [31], where superconductivity is not observed 
at zero field but emerges when a large (> 15 T) field is applied. This oc- 
curs due to the Jaccarino-Peter effect [32] which can be described using the 
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Fig. 1.2. A schematic pressure-temperature phase diagram of «-(ET)2-X\ The 
first order Mott transition is shown as a thick line, while second order transi- 
tions are thin lines and crossovers are dashed lines. The effect 'chemical' pressure 
for different anions, X, is indicated by the arrows on the abscissa. P = corre- 
sponds to ambient pressure for X =Cu[N(CN)2]Cl (abbreviated as Hg-Cl in the 
diagram). Hs-NCS indicates the effective chemical pressure for X =Cu(NCS)2 (rel- 
ative to that for X =Cu[N(CN)2]Cl) and Hs-Br indicates the effective chemical 
pressure for X =Cu[N(CN)2]Br. Dg-Br indicates the effective chemical pressure of 
X =Cu[N(CN) 2 ]Br with the ET molecule fully deuterated. (Modified from [30].) 

Hubbard-Kondo model [33] , in which a Kondo term (to describe the magnetic 
anions) is added to the Hubbard model described in Sect. 1.2.2. For a recent 
review of these materials see [34]. 

1.2.4 Multiband charge transfer salts - the a phase 

The a phase salts have a complex band structure in which many bands cross 
the Fermi level [35]. Thus they cannot be described in the one-band frame- 
work which is relevant to the other organic charge transfer salts. These salts 
have many interesting properties such as charge ordering, an unconventional 
metallic state and superconductivity [15]. However, space will not permit us 
to discuss them at length here. 

1.2.5 Quasi one-dimensional charge transfer salts 

Finally some organic charge transfer salts, e.g., the Bechgaard and Fabre salt 
families, are significantly anisotropic in all three directions and are therefore 
often described as quasi-one-dimensional conductors. Typically the hopping 
integrals are an order of magnitude larger in intrachain direction than in 
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the interchain direction, while the interchain hopping integrals are an order of 
magnitude larger than those intcrplane [36] . These materials show many inter- 
esting behaviours such as spin and charge density waves, superconductivity, 
possible Luttingcr liquid phases, Mott insulating states, spin-Pcicrls states, 
antifcrromagnctism and Fermi liquid behaviour [37]. A key question is which 
of the phenomena arise from the quasi-one-dimensionality of the system and 
for which properties higher dimensional models are required to understand 
the behaviour. 

For the remainder of this review we will discuss the half-filled organic 
charge transfer salts introduced in section 1.2.2 unless it is explicitly stated 
otherwise. 



1.3 The metallic state 

The metallic state is of fundamental importance. Firstly the organic charge 
transfer salts exhibit an unconventional metallic state which is of great interest 
in its own right. However, many states of matter occur as instabilities in 
the metallic state: think, for example, of the Stoner and Cooper instabilities. 
Thus, for example, the difficulties in understanding the metallic state of the 
cuprates have greatly compounded the difficulties in explaining the origin of 
high temperature superconductivity. 

Many features of the metallic states of the organic charge transfer salts 
seem to agree remarkably well with the predictions of dynamical mean field 
theory (DMFT). In particular, a number of experiments show the features 
predicted to occur in the crossover from a Fermi liquid to a 'bad metal' de- 
scribed by DMFT. However, there remain features of the metallic state that 
arc not well understood. For example, NMR experiments show features con- 
sistent with a pesudogap suggesting similarities to the underdoped cuprates 
(see section 1.8.2). We will only give a very brief discussion of DMFT itself so 
as to focus on the results and the comparison with experiment. However, both 
a brief, accessible introduction to DMFT [38] and a more complete technical 
review [39] are available in the literature. 

1.3.1 The bad metal and dynamical mean field theory 

The low temperature metallic states of the organic charge transfer salts are 
exceptionally pure Fermi liquids. This is indicated most clearly by the observa- 
tion of quantum oscillations and angle-resolved magneto-resistance oscillations 
(AMRO) [40,41] which demonstrate that the electron transport is coherent (at 
least within the planes) and that the quasiparticle lifetime is extremely long. 
At low temperature the resistivity varies like p(T) = po + AT 2 [14, 15,42,43], 
which is the behaviour expected of a Fermi liquid when electron-electron in- 
teractions are the dominant scattering mechanism [44,45]. 
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However, the temperature dependence of the resistivity (shown in figure 
1.3) is very different from what is expected for simple metals [where p(T) 
monotonically decreases as T is lowered]. Instead at high temperatures the 
resistivity increases as the temperature is lowered, reaching a broad max- 
imum, and the quadratic temperature dependence is only observed below 
about 30 K. If the mean free path is less than the lattice spacing it is not 
meaningful to speak of a coherent wavevector which describes the electronic 
transport. This condition allows us to define the Mott minimum conductiv- 
ity [46] (also known as the Mott-Ioffe-Regel limit) which is of order 10 3 S/cm 
for k-(ET)2Cu(NCS)2- In contrast the conductivity at the peak in the re- 
sistivity is ~1 S/cm (c.f., Fig 1.3). However, below the peak the resistivity 
does decrease monotonically, thus this state is described as a 'bad metal'. Bad 
metal behaviour is also seen in the alkali doped fullerides [47], Sr 2 Ru04 [48], 
SrRu03 [49], and VO2 [50]. All of these systems are strongly correlated ma- 
terials 'near to' a Mott transition. No Drude peak is evident in the optical 
conductivity of the organic charge transfer salts above about 50 K [51]. Fur- 
ther, even at low temperatures where a 'Drude-likc' feature does appear, this 
can only be fit to the Drude form by introducing a frequency dependencies to 
the scattering rate and the effective mass [51]. On the other hand a broad, high 
frequency peak is observed at all temperatures. This is suppressed somewhat 
at low temperatures. 

The DMFT of the Hubbard model provides a description of all of this 
physics [52,54]. DMFT works by mapping the Hubbard model onto the An- 
derson model for a single magnetic impurity in a bath of conduction elec- 
trons [38,39]. This procedure is exact in the limits of infinite spatial dimensions 
or infinite lattice coordination number. For half-filled systems close to a Mott 
transition DMFT predicts a Fermi liquid at low temperatures, and a crossover 
to incoherent transport as the temperature is raised. The crossover temper- 
ature scale, T co h, is related to the destruction of Kondo screening and Fermi 
liquid behaviour with increasing temperature (above the Kondo temperature 
Tk) in the Anderson model. In the Anderson model the conduction electrons 
are strongly scattered by a (localised) magnetic impurity for T > Tk- But 
for T < Tk a quantum coherent singlet forms between the impurity and the 
conduction electrons [3]. In the DMFT of the Hubbard model for T > T co h 
the electrons are quasi-localised and the electrons on the single site treated 
exactly strongly scatter those in the bath. However, for T < T co h transport is 
coherent and the electrons only scatter one another weakly, thus Fermi liquid 
behaviour is regained. This is why, for example, the temperature dependence 
of the resistivities of the Anderson and Hubbard models arc so similar [3,39]. 
DMFT predicts that there is no Drude peak in the bad metal phase and that 
most of the spectral weight is contained in high energy features. This is be- 
cause much of the spectral weight is transferred to the 'Hubbard bands' that 
will emerge in the Mott insulating state [38,39]. This prediction is clearly 
in good agrement with the optical conductivity measurements [51] described 
above. 
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Fig. 1.3. Non-Fermi liquid behaviour in the resistivity and thermopower of half- 
filled quasi-two-dimensional organic charge transfer salts. The circular data points in 
the left panel show the non-monotonic temperature dependence of the resistivity, p, 
in «>(ET)2Cu[N(CN)2]Cl (to be contrasted with the monotonic behaviour predicted 
by band theory). At low pressure a broad maximum is observed in the resistivity 
as is predicted by dynamical mean field theory (DMFT) (shown as rhomboidal 
data points). As the pressure is increased the broad maximum is suppressed and 
eventually becomes entirely absent. DMFT reproduces these trends if pressure is 
taken to control the relative interaction strength (P ~ W/U [28]). This data is taken 
from [52]. The right panel shows the temperature dependence of the thermopower, 
S, in the highly conducting plane of K-(ET) 2 Cu[N(CN) 2 ]Br (taken from [53]). The 
details of the non-interacting band structure dictate that the thermopower is positive 
when the heat current is parallel to the a-axis and negative when the heat current 
is parallel to the c-axis. However, along both axes the thermopower has a broad 
extremum at around the same temperature as the maximum in the resistivity occurs 
in the left panel. These extrema in the thermopower are predicted by DMFT [54] 
with U ~ W, whereas one expects the thermopower to be monotonic in a system 
with weakly interacting quasiparticles. Thus these two experiments, which do not 
have a natural explanation in a Fermi liquid picture are explained by the crossover 
from a 'bad metal' at high temperatures to a Fermi liquid at lower temperatures 
described by DMFT. 

The success of DMFT in describing the transport properties and the phase 
diagram of many organic charge transfer salts down to temperatures of about 
10 K (where, for example, superconductivity becomes important) has been 
rather puzzling, until recently [11], given that these materials are quasi-two- 
dimensional and that DMFT is only expected to be a good approximation in 
the limit of high spatial dimension or high co-ordination number. However, the 
applicability of DMFT to low-dimensional systems with large frustration is 
consistent with the fact that for frustrated magnetic models a Curie- Weiss law 
holds down to a much lower temperature than for unfrustrated models [55,56]. 
Deviations from Curie- Weiss behaviour result from spatially dependent cor- 
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relations. Hence, we expect that a DMFT treatment of the Hubbard model 
on frustrated lattices will be a good approximation down to much lower tem- 
peratures than it is for unfrustratcd models. 

However, DMFT is not simply a theory for describing the conductivity: 
good agrcmcnt with experiment is found for many other properties. For exam- 
ple, measurements of the linear coefficient of the specific heat, 7, of K-(ET)a- 
Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 [57] show that the effective mass, to*, 
is several (3-5) times that predicted by band theory [58]. Similar effective 
masses are found in quantum oscillation experiments [40,41]. DMFT predicts 
this large mass enhancement. This ia a result of the strong electronic correla- 
tions. 

Experimentally, the thermopower, shown in figure 1.3, in the highly con- 
ducting (ac) plane of K-(ET) 2 Cu[N(CN) 2 ]Br is positive (negative) in the a 
(c) direction. 2 The thermopower, S, shows a broad maximum (minimum) at 
around 100 K (60 K) [53], i.e., at about the same temperature the resistivity 
shows a maximum (figure 1.3). The thermopower predicted for a weakly in- 
teracting metallic state with the band structure of K-(ET) 2 Cu[N(CN) 2 ]Br is 
~5 times smaller than that observed experimentally [25,53]. DMFT predicts 
an extremum in the thermopower at T ~ T co h for large U /W [54] and the 
large effective mass predicted by DMFT gives rise to the large thermopower 
as S/T ex to*, in good agrement with experiment. 

A strong decrease in the ultrasonic velocity is observed at temperatures 
-40 K in K-(ET) 2 Cu[N(CN) 2 ]Br [59], k-(ET) 2 Cu(NCS) 2 [59] and fs-(ET) 2 - 
Cu[N(CN) 2 ]Cl (under pressure [60]). The pressure (and chemical pressure) 
dependence of the temperature at which these anomalies are observed show 
that they occur at the same temperature as the other anomalies that, we have 
argued above, correspond to T co h- Does this then suggest that phonons play an 
important role in this physics? DMFT shows that this is not the case. DMFT 
studies of the Hubbard-Holstein model [61,62] in the limit where the electron- 
electron interactions are much stronger than the electron-phonon interactions 
(U 3> g) predict phonon anomalies at T co h- However, these anomalies are 
parasitic: the changes in the behaviour of the electrons cause the change in the 
behaviour of the lattice vibrations because of the electron-phonon coupling, 
and the phonons do not play a significant role in driving the crossover from 
the Fermi liquid to the 'bad metal'. 

The Sommerfeld- Wilson ratio is defined as R\y = [x{T — 0)70]/ (7X0), 
where x(T) is the magnetic susceptibility, xo is the zero temperature magnetic 
susceptibility of the non-interacting electron gas and 70 is the linear coefficient 
of the heat capacity for the non-interacting electron gas. Clearly for the non- 
interacting system R\y = 1, but the Kondo model predicts Rw = 2. Therefore, 
as DMFT exploits the deep connections between Kondo physics and the Mott- 



2 The sign difference results from the details of the band structure of k-(ET)2- 
Cu[N(CN)2]Br and corresponds to hole (electron) like conduction along the a (c) 
directions [53]. 
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Hubbard transition DMFT predicts R\y > 1 in the organics. Experimentally 
R w is 1.5 ± 0.2 in both K-(ET) 2 Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 [63]. 

The Kadowaki- Woods ratio is Rkw = A/^ 2 , where A is the quadratic 
coefficient of the resistivity. In many strongly correlated materials Rkw ~ 
10.5 u,Q cm K 2 mol 2 J -2 . The Kadowaki- Woods ratio is significantly larger 
than this in the half-filled layered organic materials, even when differences in 
the unit cell volume are allowed for [64,65]. This has led to suggestions that 
the quadratic temperature dependence of the resistivity may not result from 
strong electronic correlations, but from phonons [14,66]. However, one must be 
careful to allow for the fact that the only reliable resistivity measurements in 
the organic charge transfer salts are perpendicular to the highly conducting 
plane. When this is allowed for the observed Kadowaki- Woods ratio is the 
expected order of magnitude [65] . 

Collectively the experiments described above (and those studying the Mott 
transition which we review in section 1.4) show that a weakly interacting, 
Fermi liquid, description is not sufficient to explain the full temperature de- 
pendence of the thermodynamic and transport properties of the half-filled 
organic charge transfer salts. However, DMFT, which includes the effects of 
strong electron-electron interactions, and reproduces Fermi liquid theory be- 
low T co h, can provide both a quantitative [52] and qualitative [54] description 
of the full temperature dependence. However, DMFT is a purely local theory. 
Therefore, any properties with a significant k-dependence will not be properly 
described by DMFT. We will discuss some such possible features in section 
1.8.2. 



1.4 The Mott transition 

DMFT also provides a description of the, first order, Mott metal-insulator 
transition [38,39]. For U = the density of states (DOS) calculated from 
DMFT is, as one should expect, that of the tight-binding model. For a weakly 
interacting system the DOS will be only weakly renormalised by the interac- 
tions. However, as U is turned on the system becomes strongly interacting and 
a peak in the DOS will emerge at the Fermi energy. This peak is associated 
with quasiparticles. The total spectral weight associated with the quasipar- 
ticles is Z = mb/m* < 1, where nib is the band mass of the electron, hence 
Z is called the quasiparticlc weight. The remainder of the spectral weight is 
transferred to two broad bands centred on ±[//2 of width ~ W, known as the 
Hubbard bands. The Hubbard bands correspond to quasi-localised states with 
the lower Hubbard band (centred at —U/2) corresponding to singly occupied 
sites and the upper Hubbard band (centred at +U /2) corresponding to doubly 
occupied sites. As U is further increased the quasiparticlc peak narrows (Z 
decreases) and more spectral weight is transferred to Hubbard bands. When 
U is increased above the critical U ~ W for the Mott transition the quasi- 
particle peak vanishes and all of the spectral weight resides in the Hubbard 
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bands. There is now no density of states at the Fermi energy and so we have 
an insulator. A number of studies have found that this transition is first order 
within the DMFT framework [39,67]. 

It is conceptually useful to compare the Mott transition in the organic 
charge transfer salts (figure 1.2) and that predicted for the Hubbard model 
on the anisotropic triangular lattice by DMFT with the liquid-gas transition. 
The difference between the liquid and the gas is their densities. One may move 
from the gas to the liquid in cither of two ways: either directly, by increasing 
the pressure through a first order phase transition where the density changes 
discontinuously; or by first increasing the temperature above the critical tem- 
perature and then passing around the critical point. If one passes around the 
critical point no phase transition is observed and the density varies continu- 
ously through the fluid phase. The Mott transition behaves in much the same 
way. In the Mott insulator the electrons are localised and hence the conductiv- 
ity is poor whereas the metal has a large conductivity. We may either drive the 
Mott transition by increasing the pressure and passing directly through the 
first order transition where the degree of localisation changes discontinuously, 
or else we may first increase the temperature above the critical temperature 
where the line of first order transitions ends (see figure 1.2) and thus pass 
continuously from the insulating phase to the metallic phase. In this picture 
the 'bad-metal' is somewhat analogous to the fluid. As we move through the 
'bad-metal' the conductivity changes continuously as we pass from the lo- 
calised behaviour of the insulator to the coherent excitations characteristic 
of the Fermi liquid. Although one should note that the analogy is not exact 
as the 'bad metal' regime extends slightly below the critical point (c.f. figure 
1.2). 

1.4.1 The critical point 

We now turn our discussion to the critical point itself. Unsurprisingly, DMFT 
predicts mean-field critical exponents (see table 1.1). Limclettc et al. [68] 
measured the critical exponents of the Mott transition in V2O3 [68]. 3 They 

3 These exponents can be measured by studying the conductivity near the Mott 
transition. The conductivity gives access to the critical behaviour as it is related to 
the thermodynamic ground state via linear response theory. The critical exponents 
can be determined from the conductivity, cr, as: (a — a c ) ~ (T c — T) 13 , (da /<IP)t ~ 
\T — T c |~ 7 and (<r — a c ) ~ (P — P c ) 1//<s , where a c is the conductivity at the critical 
point, T is the temperature, T c is the critical temperature, P is the pressure and 
P c is the critical pressure. To interpret the order parameter exponent on the critical 
isotherm, 8, one should recall that the order parameter of the Mott transition is 
the half-width of the coexistence curve, which vanishes at the critical point. This 
is entirely analogous to the liquid-gas transition [69]. It is often stated that the 
order parameter of the liquid-gas transition is the difference in the densities of the 
liquid and the gas. This statement needs to be interpreted carefully. Strictly the 
density difference between the liquid and the gas phase when they are in coexistence 
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Table 1.1. Comparison of the critical exponents observed for some common univer- 
sality classes with those seen at the critical point of the Mott transition in k-(ET)2- 
Cu[N(CN)2]C1. The symbol X — > Y indicates that the measured exponent crosses 
over from X in a small critical region near the critical point, to Y in a larger region 
further from the critical point. 



observed three dimensional (3D) Ising exponents (the same as are seen in the 
liquid-gas transition), but only in a extremely narrow critical region around 
the critical point (P/P c , T/T c < 1CP 2 ). Limclettc et al. observed mean-field 
critical exponents except in this extremely small critical region near the crit- 
ical point which, they argued, implies that there are bound pairs of doubly 
occupied and vacant sites near the Mott transition. These bound states can 
exist over large length scales (of order nm) and are argued to be the root 
cause of the success of DMFT. This is analogous to the success of the (mean 
field) BCS theory of superconductors which results from the large coherence 
length (and hence small critical region) in conventional superconductors. 

However, recently Kagawa et al. [70] have found that the critical exponents 
for the critical end point of the Mott transition in «-(ET)2Cu[N(CN)2]Cl do 
not correspond to either the Ising or mean field universality classes. Indeed the 
critical exponents do not belong to any universality class that had previously 
been observed (see table 1.1). Crucially, in spite of being very different from 
those for the Ising, XY or Heisenberg universality classes, the observed expo- 
nents obey the standard scaling relation (5 — l)f3 = 7, within experimental 
error. 

Remarkably, Imada [71] had predicted these exponents from phenomcno- 
logical theories of the Mott transition. Misawa et al. [72] have now shown that 
these critical exponents may be derived within the Hartree-Fock approxima- 
tion from the Hubbard model on the anisotropic triangular lattice. Within 
the Hartree-Fock approximation the metal becomes magnetically ordered at 



is a good order parameter. This is, of course, proportional to the half-width of the 
coexistence curve because of the Maxwell equal area construction. However, the 
density difference when the liquid and gas when they are not in coexistence is not 
a good order parameter. These comments are equally valid for the difference in 
conductivities between the metal and insulator in the Mott transition. 
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relatively small U/W. But the metal insulator transition does not occur until 
U ~ W. Thus there is no symmetry breaking at the metal- insulator tran- 
sition within the Hartree-Fock approximation. Misawa et al. find that for 
0.056 m t' cl /t < t'/t < t' c2 /t ss 0.365 a line of first order phase transitions 
ends at a finite temperature critical point. At t' cl and t' c2 this critical point is 
driven to T = and for t'/t < t' cl /t or t'/t > t' c2 /t a quantum critical point is 
observed, t' = t' cl and t' = t' c2 are therefore termed marginal quantum critical 
points (MQCPs). At the MQPCs Misawa et al. find that the critical exponents 
are (3 = 1, 7 = 1 and 5 = 2, which satisfy the scaling relation (5 — l)j3 = 7, 
and agree with those measured by Kagawa et al. [70] in k-(ET) 2 Cu[N(CN) 2 ]C1 
within experimental error. However, one should recall that these exponents 
are only strictly valid at T = and the critical point in k-(ET) 2 Cu[N(CN) 2 ]C1 
occurs at 39.7 K - where we might expect to regain the standard Ising critical 
exponents. However, Misawa et al. stress that the true critical exponents will 
only apply within a small critical region very close to the phase transition and 
that beyond that region the observed exponents will crossover to those char- 
acteristic of the MQCP. Misawa et al. estimate that experimental resolution 
of Kagawa et al. is two orders of magnitude worse than that required to see 
the Ising critical behaviour. Misawa et al. stress that the fact that magnetic 
order exists on both the insulating and metallic sides of the phase transition 
is vital for their theory 4 We should note that that in k-(ET) 2 Cu[N(CN) 2 ]C1 
antiferromagnetic ordering does not occur over the entirety of the Mott in- 
sulating phase, but only at low temperatures (see figure 1.2). Therefore, the 
Mott transition in k-(ET) 2 Cu[N(CN) 2 ]C1 is not accompanied by any symme- 
try breaking in the vicinity of the critical point. Importantly, the Ginzburg 
criterion shows that the upper critical dimension is two and so the mean field 
Hartree-Fock treatment, which neglects all fluctuations, appears likely to be 
correct in cases where no symmetry is broken at the critical point. It is im- 
portant to note the role that frustration plays in this scenario as for weakly 
frustrated lattices (t'/t < t' cl /t or t'/t > t' c2 /t) these unconventional critical 
exponents are not predicted. 



1.5 The superconducting state 

The organic charge transfer salts are often referred to as the organic super- 
conductors. However, this is somewhat misleading both because of the range 
of other phenomena observed in the organic charge transfer salts and because 
of the range of other materials which fit the description 'organic superconduc- 
tors'. Other organic superconductors include the alkali doped fullcridcs [47], 
intercalated graphite [73] and ion implanted polymers [74]. However, below we 



4 Note that Imada's [71] phenomenological theories which predict the same ex- 
ponents do so in the case of the non-magnetic-insulator to non-magnetic-metal 
transition. 
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discuss only superconductivity in the quasi-two-dimensional half-filled organic 
charge transfer salts. 

The phrase 'unconventional superconductivity' has several meanings. These 
meanings are best illustrated by comparison to elemental superconductors, 
which are described by BCS theory [4] (or more accurately by Eliashbcrg the- 
ory [75] or the density functional theory of superconductivity [76]). In these 
theories the electrons form Cooper pairs due to an effective attractive interac- 
tion between electrons mediated by the exchange of phonons. The BCS state 
breaks gauge symmetry but no other symmetry of the crystal. Thus the phrase 
'unconventional superconductivity' can be applied to any superconductor that 
does not satisfy any one of these three requirements (BCS, phononic pairing 
mechanism and no additional symmetry breaking) for conventional supercon- 
ductivity. In this section we review the experimental evidence which shows 
that the superconductivity in the half-filled organic charge transfer salts is 
unconventional in all three of the senses described above. We delay a detailed 
discussion of the theory until section 1.6. 

1.5.1 Pairing symmetry (evidence for additional symmetry 
breaking) 

In the original formulation of BCS theory [4] one assumes that the ef- 
fective attractive pairwise interaction responsible for superconductivity, V, 
is spatially uniform. This means that the superconducting order parame- 
ter A is also isotropic as A = X)k ^( c kT c -kl)- This assumption is some- 
what unphysical, but it is straightforward to generalise BCS theory to al- 
low for momentum dependent effective interactions, Vk [77]. In this case we 
should also allow the order parameter to develop a k-dependence by defining 
A^ = J^ k , Vk-k'(ck'T c -k'i)- It is now natural to ask what symmetry A^ has. 
Note that although we have introduced these ideas in the context of BCS the- 
ory the symmetry arguments that follow arc based basic truths about quantum 
mechanics and do not depend on the details of the microscopic theory of super- 
conductivity. Thus symmetry analyses are a powerful framework in which to 
understand the phenomenology of unconventional superconductors [7,24,78]. 

A brief introduction to the symmetry of pairing states 

A fundamental theorem of quantum mechanics is that the eigenstates must 
transform according to an irreducible representation of the symmetry group of 
the Hamiltonian [79] . 5 This is true regardless of the complexity of the quantum 

5 Here we neglect the possibility of accidental degeneracies as they complicate 
the discussion considerably, but do not really change the physics as the accidentally 
degenerate states can always be decomposed into a set of eigenstates that trans- 
form under the operations of the group in the same way as a set of the irreducible 
representations of the symmetry group of the Hamiltonian. 
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many-body Hamiltonian. Therefore, one expects on very general grounds that 
the superconducting order parameter (which has the same symmetry as the 
wavefunction describing the superconducting state 6 ) will have the symmetry 
of a particular irreducible representation of the group describing the symmetry 
of the normal state (and hence the Hamiltonian at T c ). This argument can be 
made rigourous exactly at the critical temperature, and may be expected to 
hold below T c provided there are no further phase transitions. 

In elemental superconductors it is found that, although Z\k is not com- 
pletely isotropic [80] the order parameter does have the same symmetry as 
the crystal. The first material found where this did not appear to be the case 
was in superfluid 3 He. Clearly there is no crystal in liquid 3 He and so the 
normal state has the full symmetry of free space. Therefore, the expectation 
that the order parameter will transform like a particular representation of the 
group describing the symmetry of the normal state can be recast as the claim 
that if we expand the order parameter in terms of the spherical harmonics, 

OO I 

Ac = ^ E VlmYim(k), (1.7) 

1=0 m=-l 

then we will find that we only require r\i m to be finite for one particular value 
of I. Thus we find that 

l 

= Vi m Y lm {k). (1.8) 

m=—l 

It is natural to refer to superconductors in which the non-zero r\i m arc I = 
as s-wave, superconductors in which the finite i]i m are I = 1 as p-wave, 
superconductors in which the finite rji m are / = 2 as d-wave, and so on, by 
analogy with atomic physics. It has been established that 3 He is a p-wave 
superconductor [81]. 

Historically, the discovery that superfluid 3 He has a lower symmetry than 
normal 3 He proceeded any analogous discoveries in superconductors by some 
time. However, there are now several materials in which the superconducting 
state is believed to have a lower symmetry than the normal state. These in- 
clude, Sr2Ru04, the cuprates, several heavy fermion materials and, we will 
argue below the layered organic charge transfer salts. One often hears state- 
ments such as "the cuprates are <i-wave superconductors". It is not imme- 
diately clear what this means as the normal state of the cuprates does not 
have the full symmetry of free space (due to the crystal lattice) and so the 
spherical harmonics are not an appropriate basis in which to expand the order 
parameter. 

6 Assuming no other phase transition accompanies the superconducting transi- 
tion. 
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In a crystal the natural basis is that of the irreducible representations of the 
point group of the crystal. 7 (Table 1.2 describes the point group which is 
that of most cuprates and a number of other unconventional superconductors 
and tables 1.3-1.5 give three examples of point groups relevant to the organic 
charge transfer salts.) In this case the order parameter may be written as 



where ip[ (k) are the basis functions of the r th irreducible representation and 
d r is the dimension of r. 

The above discussion is only valid for singlet superconductors. Triplet su- 
perconductivity does not greatly complicate this analysis, but one does have 
to include a sum over the three spin projections of the S = 1 Cooper pairs 
and generalise the order parameter slightly (sec [7], [78] or [82] for a careful 
discussion of the group theoretic classification of triplet states). 

Thus when the statement is made that "the cuprates are d x i_ y i super- 
conductors" , what is meant is that the superconducting state transforms like 
the Big representation of the group (c.f. table 1.2). However, all possi- 
ble basis functions of the B\ g representation of vanish along the lines 
k x = ±k y . Thus the nodes of the order parameter (which are often all one is 
concerned with, see below and section 1.5.2) are the same as for the d x 2_ y 2 
spherical harmonic. Thus in analogy with the case of 3 He, people often re- 
fer to 's-wave', 'p-wave' or 'd-wave' superconductors. One should always bare 
in mind what this naming conventional really means, as confusion can arise 
(and repeatedly has arisen) when the terms are used carelessly. However, this 
terminology is all but universal and therefore we will make use of it, but we 
will use inverted commas to remind the reader that an analogy is being drawn 
and that we really mean a particular irreducible representation of the point 
group of the crystal in question. 

Singlet or triplet? 

As all of the half-filled organic charge transfer salts have inversion symmetry 
singlet and triplet states are distinct [78]. Thus the first question we should 
ask is whether the superconducting states of these materials are singlet or 
triplet. Measurements of the 13 C NMR Knight shift [83-85] in k-(ET) 2 Cu- 
[N(CN)2]Br, with the magnetic field, H, parallel to the conducting planes, 
show that as T — > so does the Knight shift. This single experiment does 
not actually rule out triplet pairing, although it does make triplet pairing 
extremely unlikely. This experiment is compatible with a triplet state in which 
d(k) x H = where d(k) is the usual Balian-Werthamer order parameter 

7 We assume throughout that the superconducting state does not break transla- 
tional symmetry. 
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Table 1.2. The character table, symmetry required nodes in the superconducting 
energy gap and some basis functions of the irreducible representations of the point 
group Dih- This is the point group of many unconventional superconductors with 
tetragonal crystals, including many of the cuprates and Sr2Ru04- We assume that 
the x, y and z axes are, respectively, parallel to the a, b and c axes. The functions 
lk, X k , Yk and Z k may be any functions which transform, respectively, as 1, k x , 
k y and k z under the operations of the group and satisfy translational symmetry. 
The operations of the group are the identity (E), rotation by n and tv/2 about 
the c axis (respectively Cf and C4), rotation by n about the a and b axes (which 
have the same characters, which we therefore label C^ 6 ), rotation by ir about the 
either of the diagonal of the a-b plane (which have the same characters and are 
labelled Cf), inversion (i: inversion symmetry takes k to — k), an improper by 7r/4 
about the c axes {S%: an improper is a rotation followed by a reflection through the 
plane perpendicular to the axis of rotation), reflection through the ab plane (a ab ) 
reflection through the ac and be planes (which have the same characters, which we 
therefore label a~^) and reflection through the planes specified by (x + y)z = 
and (x — y)z — (which have the same characters and are labelled a dc ). A brief 
explanation of characters is given in the caption to table 1.3. 

for triplet superconductivity [82,86]. However, Zuo et al. [87] measured the 
critical field as a function of temperature with H parallel to the conducting 
planes. In this configuration no orbital currents flow so the critical field is due 
to Clogston-Chandrasckhar (or Pauli) limit [88-90]. There is no Clogston- 
Chandrasekhar limit for H _L c for triplet states compatible with measured 
Knight shift. Thus, for such states there would be no critical field with H||6 (in 
fact for such states one would increase T c by applying a field parallel to the b- 
axis [90]). Experimentally [91] it is found that superconductivity is destroyed 
by a magnetic field parallel to the b-axis. Therefore, only when considered 
together do the three experiments discussed above [83,87,91] strictly rule out 
triplet pairing. Further evidence for Clogston-Chandrasckhar limiting, and 
hence singlet pairing, comes from the observation that the in plane upper 
critical field is independent of the field direction [92]. Given the anisotropic 
nature of the Fermi surface of K-(ET)2Cu[N(CN)2]Br it is extremely unlikely 
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Table 1.3. The character table, symmetry required nodes in the superconducting 
energy gap and some basis functions of the irreducible representations of the point 
group Z?2Ji- This is the symmetry of the orthorhombic organic superconductors such 
as K-(ET)2Cu[N(CN)2]Br. We argue in section 1.5 that the experimental evidence 
(see Figs. 1.4 and 1.5) shows that the superconducting order parameters of fc-(ET)2- 
Cu[N(CN)2]Br and other orthorhombic organic charge transfer salts transform like 
the Big representation of Dm- Thus the superconducting state may be said to be 
l d x 2_ y 2 \ In the colloquial names column we include d xy , parenthetically, as a name 
for the Aig representation. This is not intended to encourage the use of this name 
but merely to stress that order parameters transforming like the Big representation 
of D^h (which have been invoked to explain some experimental data; see, e.g. [93]) 
transform according to the Ai g representation of D2h and are therefore are not sym- 
metry distinct from 's-wave' order parameters. Note that in K-(ET)2Cu[N(CN)2]Br 
the highly conducting plane is the a-c plane, and that the x and y axes are taken 
to lie along the same directions as the t hopping integrals (c.f., figure 1.1) therefore 
k a — kx + ky and k c = k x — k y . The functions lk, ^k, Yk, Zk, Ak, Bk and C k may be 
any functions which transform, respectively, as 1, k x , k y , k z , k a , kb and k c under the 
operations of the group and satisfy translational symmetry. The operations of the 
group are the identity {E), rotation by n about the a, b and c axes (respectively , 
C\ and C 2 a ), inversion (i) and reflection through the ab, ac and be planes (respec- 
tively a ab , a ac and a bc ). The character is the trace of any matrix that can represent 
the operation in that irreducible representation (see [79] for a detailed discussion). 
However, for the one-dimensional representations relevant here the character has 
more physical interpretation: the character is the sign introduced in the order pa- 
rameter (and wavefunction) by that operation. For example, any order parameter 
transforming according to the A\ g representation will by unchanged by any of the 
operations of the group. On the other hand if the order parameter transforms ac- 
cording to the B2g representation (as does, for example, At ~ cosk x — cosk y ) then 
the order parameter changes sign under rotation by n about the a (x + y) or c (x — y) 
axes and reflection through the be (x = y) and ab (x — —y) planes. 

that orbital mechanisms for the destruction of superconductivity would be so 
isotropic. 



1 Strong electronic correlations in organic charge transfer salts 21 
Spatial symmetry ('s-wave' versus 'd-wave') 

The organic charge transfer salts form orthorhombic, monoclinic and triclinic 
crystals. This can lead to important differences between their superconduct- 
ing states [24,94]. Orthorhombic crystals such as K-(ET) 2 Cu[N(CN) 2 ]Br (in 
which the highly conducting plane is the a-c plane) are described by the D2h 
point group. There are four irreducible representations of D^k which corre- 
spond to singlet superconductivity (see table 1.3). Given the layered struc- 
ture of the crystals superconductivity transforming as either the B\ g or B^, g 
irreducible representations is unlikely [24,78]. Therefore, our task is to differ- 
entiate between superconducting order parameters transforming like the Ay g 
representation (which is often referred to as 's-wave' superconductivity) and 
superconducting states that transform according to the Bi g representation 
(or l d x 2_ y 2^ superconductivity). 

There are even less symmetry distinct superconducting states in mono- 
clinic crystals such as k-(ET)2Cu(NCS) 2 (for which the highly conducting 
plane is the b-c plane) which are described by the Cih point group. There 
are two irreducible representations of Cih which correspond to singlet su- 
perconductivity (see table 1.4). Therefore, the singlet order parameter must 
transform like either the A g representation (which we may refer to as 's-wave' 
superconductivity) or the B g representation (which wc may refer to as 'd-wave' 
superconductivity 8 ) . 

One way to differentiate between 's-wave' and 'd-wave' states 9 is to mea- 
sure the low temperature behaviour of thermodynamic variables. For an 's- 
wave' state the superconducting gap is finite at every point on the Fermi 
surface. As excitations may only take place at energies within about fc^T of 
the Fermi energy thermodynamic properties arc exponentially activated at 
low temperatures. For example, the heat capacity C v oc exp(— Ao/ksT) for 
T -C T c , where Aq is the minimum value of the magnitude of the supercon- 
ducting gap at the Fermi surface at T = 0. In contrast for a l d x 2_ y 2' state 
symmetry requires that the superconducting order parameter vanishes along 
(at) four lines (points) on a three (two) dimensional Fermi surface. These 
lines (points) are known as nodes. At the nodes there are excitations with 
arbitrarily low energies. Thus although density of states, D(e), does vanish at 
the Fermi energy, ep, in a 'd-wave' superconductor, the DOS grows linearly 
as we move away from the Fermi energy, D(e) oc \e — ep\ for \e — ejr | <C \A\. 
As, for example, C v /T ~ D{t) it follows that at low temperatures C v oc T 2 in 
a 'd-wave' superconductor. 

Measurements of the temperature dependence of the heat capacity [13], 
penetration depth [95] and nuclear spin lattice relaxation rate [83-85] have 

8 Although the only symmetry required node for states transforming according 
to the B g representation lies along the c-axis [24] and so such states will not, in 
general, be the same as the i d x i_ y i' state described above 

9 From hereon we use the phrase 'd-wave' to refer to both the B g representation 
of Cnh and the B2 g representation of D^h unless we explicitly state otherwise. 



22 B. J. Powell and Ross H. McKenzie 



Irrep 




i 




Required 
nodes 


Example basis functions 


Colloquial names 


A g 


1 1 
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1 


none 


Ik, A k , B£, C^, X k Y k 






B a 


1 -1 


1 


-1 


line 




Y 2 

f k 


d 


A u 


1 1 


-1 


-1 


none 


A k , Bk, X k + Yk 




Px + y 


B u 


1 -1 


-1 


1 


line 


Ck, Xk — Vk 




Px — y 



Table 1.4. The character table, symmetry required nodes and some basis functions 
of the even parity irreducible representations of the point group C2J1. This represents 
the symmetry of the monoclinic organic superconductors such as k-(ET)2Cu(NCS)2- 
We argue in section 1.5 that the experimental evidence (see Figs. 1.4 and 1.5) shows 
that the superconducting order parameters of k>(ET)2Cu(NCS)2 and other mon- 
oclinic organic charge transfer salts transform like the B g representation of C2h- 
Thus the superconducting state may be said to be 'd-wave'. In the colloquial names 
column we include d xy , parenthetically, as a name for the Ai g representation. This 
is not intended to encourage the use of this name but merely to stress that order 
parameters transforming like the Bi g representation of (which have been in- 
voked to explain some experimental data; see, e.g. [93]) transform according to the 
Ai g representation of Cih are therefore are not symmetry distinct from 's-wave' 
order parameters. Note that in k-(ET)2Cu(NCS)2 the highly conducting plane is 
the b-c plane, and that the x and y axes are taken to lie along the same directions 
as the t hopping integrals (c.f., figure 1.1) therefore kb = k x + k y and k c = k x — k y . 
The functions lk, X k , Yk, A k , B k and Ck may be any functions which trans- 
form, respectively, as 1, k x , k y , k z , k a , kt and k c under the operations of the group 
and satisfy translational symmetry. A brief explanation of characters is given in the 
caption to table 1.3. 



produced apparently rather contradictory results; with some groups arguing 
that their results provide evidence for fully gapped ',s-wave' pairing and other 
groups arguing that their data favours states with nodes such as 'd-wave' su- 
perconductivity. (We recently reviewed this data in rather more detail in [96].) 
However, it is important to note that many of these experiments were not 
performed in the T -C T c limit and therefore do not strongly differentiate 
between exponentially activated and power law behaviours. The only exper- 
iment, which we are aware of, reporting data below about 20% of T c is the 
penetration depth, A, study of Carrington et al. [95] which reports data down 
to about 1% of T c . They found a power law behaviour, but with a peculiar 
A(T) — A(0) oc T 3 / 2 behaviour 10 (see figure 1.5). The correct interpretation of 
this temperature dependence is not clear at the current time, except to note 
that their data is quite inconsistent with an 's-wave' order parameter. 

In the cuprates phase sensitive probes, especially tunnelling [98] and scan- 
ning tunnelling microscopy (STM) [99] experiments have been particularly 
important for determining the pairing symmetry. However, such experiments 



10 One expects A(T) — A(0) oc T for line (point) nodes and a three- (two-) dimen- 
sional Fermi surface and A(T) — A(0) oc T 2 for point nodes on a three-dimensional 
Fermi surface [7]. 
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Fig. 1.4. Experimental evidence for unconventional superconductivity in half- filled 
quasi-two-dimensional organic charge transfer salts from measurements of suppres- 
sion of T c by disorder (taken from [43] ) . It shows the variation of the critical temper- 
ature, T c , and residual resistivity, po, of ft-(ET)2Cu(NCS)2 as disorder is introduced 
by irradiating the sample. The line shows the prediction of the Abrikosov-Gorkov 
theory (1.10) for a 'non-s-wave' superconducting order parameter (that transforms 
as a non-trivial representation). In contrast for an 's-wave' order parameter Ander- 
son's theorem [97] predicts that the critical temperature is not suppressed by low 
levels of disorder. Thus, the strong initial suppression of T c as po (which is inversely 
proportional to the electron-impurity scattering rate) increases is strongly sugges- 
tive of a 'non-s-wave' order parameter. However, it is not currently understood why 
T c is greater than is predicted by the Abrikosov-Gorkov theory for po > 2 J7cm. 

have not been reliably performed in the organic charge transfer salts (sec [96] 
for a critical review of the experiments that have been performed to date.) 

An important distinction between 's-wave' superconductors and 'non-s- 
wave' superconductors is the effect of non-magnetic disorder. When a quasi- 
particlc is scattered by a non-magnetic impurity its momentum is changed by 
a random amount. However, Fermi statistics dictate that quasiparticles may 
only scatter from and to states near the Fermi surface. Therefore, in a super- 
conductor impurity scattering has the effect of averaging the order parameter 
over the Fermi surface [96, 100]. For states transforming as the trivial repre- 
sentation (e.g., the A\ g representation of D^h, see table 1.3) the average of 
the order parameter over the Fermi surface will, in general, be non-zero, and 
hence T c is not suppressed [97]. However, for order parameters which trans- 
form under the operations of the point group like any representation other 
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Fig. 1.5. Experimental evidence for unconventional superconductivity in half-filled 
quasi-two-dimensional organic charge transfer salts from the temperature depen- 
dence of the penetration depth (relative to the penetration depth extrapolated to 
T = 0). The data is taken form [95]. Samples a and b are K-(ET) 2 Cu[N(CN) 2 ]Br, 
while samples c and d are k-(ET)2Cu(NCS)2 (note that the non-linear scale on the 
abscissa and that the data have been offset for clarity). The inset shows the data for 
samples a and d on a log scale. It is clear that AX ~ (T/T c ) 5 . For an 's-wave' gap 
(transforming like the identity representation) one expects AX ~ exp(— ksT/Ao) 
at low temperatures whereas one expects a power law behaviour when there are 
nodes in the gap. However, there remains an unexplained feature in this data: the 
power law is not what is expected in simple theories where one expects AX to vary 
either linearly or quadratically with temperature at low temperatures. Experiments 
such as this and that in figure 1.4 suggest that the pairing state of the organic 
charge transfer salts has a lower symmetry than the crystal. In particular we argue 
in section 1.5 that orthorhombic materials, such as K-(ET)2Cu[N(CN)2]Br, have or- 
der parameters which transform like the B-2, a representation of D-2.H and so may be 
termed 'd x 2_ y 2' superconductors. We also argue that monoclinic materials, such as 
k-(ET)2Cu(NCS)2, have order parameters which transform like the B g representa- 
tion of C2/1 and so they may be termed 'd-wave' superconductors. 

than the trivial representation, that is for all 'non-s-wave' states, the aver- 
age of the order parameter over the Fermi surface must vanish by symmetry. 
Therefore, it can be shown [100, 101] that non-magnetic disorder suppresses 
the critical temperature of unconventional superconductors according to the 
Abrikosov-Gorkov [102] formula, 
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where T c q is the superconducting critical temperature in the pure system, 
tp{x) is the digamma function and r is the quasiparticlc lifetime due to scat- 
tering from impurities. It was recently pointed out that a large number of 
experiments in the literature show that in the organic charge transfer salts 
presumably non-magnetic impurities suppress T c in just the way predicted by 
the Abrikosov-Gorkov formula [96]. This has led to new experiments specif- 
ically designed to investigate the role of disorder in organic charge transfer 
salts [43]. These experiments have, however, not conclusively resolved the 
pairing symmetry. Indeed these experiments have resulted in a new puzzle. 
Analytis et al. [43] irradiated samples of k-(ET)2Cu(NCS)2 with both protons 
and x-rays (both leading to similar results). Fig 1.4 shows the observed critical 
temperature plotted against the residual resistivity of the sample. While this 
initially follows the Abrikosov-Gorkov curve, for larger irradiation doses, the 
suppression of T c is weaker than predicted. While this behaviour seems rather 
inconsistent with a pure 's-wave' order parameter a detailed explanation of 
this phenomena is lacking. 

The vortex lattice in the Abrikosov phase can yield important information 
about the pairing symmetry. In particular for an isotropic order parameter one 
expects a triangular vortex lattice [103,104]. While, for an anisotropic order 
parameter cither a square or a triangular lattice may occur as the energy 
difference between the square and triangular vortex lattices is less than 1% 
even for an isotropic order parameter [103,104]. In particular square vortex 
lattices are found in Sr2Ru04, UPt3 and YP^CuaOy-z [104], all of which are 
believed to be unconventional superconductors. In the organic charge transfer 
salts a triangular lattice is observed [105] by muon spin relaxation (/xSR) 
experiments, which does not give a strong indication of what the pairing 
symmetry is. 

We believe that a l d x 2_ y 2^ order parameter transforming as the Bi g repre- 
sentation of D2h in the orthorhombic materials and a 'd-wave' order parameter 
transforming as the B g representation of Cih in the monoclinic crystals [24] 
is most consistent with currently available data. However, there is no clear 
'smoking gun' experiment supporting this conclusion. We eagerly await such 
an experiment and in the next section we review some of the further experi- 
ments that could be used to probe the pairing symmetry. 

1.5.2 Possible probes of the pairing symmetry 

We will now briefly discuss potential experiments that might provide an unam- 
biguous determination of the pairing symmetry. Particular attention should 
be paid to directional probes such as ultrasonic attenuation, thermal con- 
ductivity and experiments on single crystals in a magnetic field. Directional 
probes have yielded important information about the gap structure in the 
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cuprates [98], UPt 3 [106] and Sr 2 Ru0 4 [107], so we will review what has been 
learnt by such methods these materials. 

Effect of an orientated magnetic field in the conducting plane on 
the heat capacity, nuclear spin lattice relaxation rate, penetration 
depth and interlayer resistance 

In the vortex state of a superconductor with nodes in the energy gap it can be 
shown that one effect of the applied magnetic field is to introduce a Doppler 
shift to the quasiparticle energy due to circulating supercurrents [108,109]. 
Using this semiclassical theory Vekhter et al. [110] showed that in a l d x 2_ y 2' 
superconductor with a circular Fermi surface the density of states at the Fermi 
level, D(Ef, a), where a is the angle between the field and the antinodal direc- 
tion, has a four fold variation as a magnetic field is rotated around the highly 
conducting plane. D(EF,ct) is maximal when H is aligned in the antinodal 
direction and minimal for H parallel to the nodes. 

The variation in D(Ep, a) should be directly reflected in the behaviour of 
the electronic contribution to the specific heat C v (a) which is predicted [110] 
to show a four-fold variation as H is rotated around the highly conducting 
plane. Further Vekhter et al. predict that C v /T oc \/HT for the field in the 
nodal direction and that C v /T oc ^/~HT 2 for the field in the antinodal di- 
rection. Similar effects are predicted for the H(a) dependence of the nuclear 
lattice relaxation rate 1/T^T, the penetration depth A [110] and the interlayer 
resistance [111]. 

Vekhter et al. also discuss the effect of lowering the symmetry of the unit 
cell and thus introducing a small 's-wave' admixture to the order parameter. 
They show that even a very small 's-wave' component reduces the four-fold 
variation of D{Ep,a) to a two-fold variation. However, experimentally great 
care would be required to detect a two-fold variation in, for example, C v as 
any slight mismanagement of the field out of the plane would cause a similar 
two- fold variation because of the anisotropy between H C 2\\ and i? C 2j_, where 
H C 2\\ is the in plane upper critical field and H C 2j_ is the out of plane upper 
critical field [112]. 

Deguchi et al. [112] have observed four-fold variation in C v (a) in S^Ru- 
O4. However, they interpret this as evidence of an anisotropic (but nodeless) 
gap on the 7 sheet of the Fermi surface. In fact Deguchi et al. observed two 
different four-fold observations as a function of the applied magnetic field 
strength. For |H| < H c2 \\ they observe an oscillation which they attribute to 
the anisotropy of H c2 n in the basal plane of S^RuO,}. At intermediate fields 
there are no oscillations in C v (a) and for H c \ < |H| <C H c2 \\ Deguchi et al. 
observe oscillations in C v {a) that are 7r/4 out of phase with those in the high 
field region. It is these low field oscillations that appear to associated with the 
anisotropy of the superconducting gap. In the organic charge transfer salts the 
in plane upper critical field is Clogston-Chandrasekhar (or Pauli) limited, [87] 
and there is therefore no anisotropy in H C 2\\{a) so one should not expect the 
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four-fold oscillations in the high field region because of variations in H C 2\\ (a). 
However, this observation does indicate the ease with which extraneous effects 
can enter this type of experiment and therefore the care with which any such 
data needs to be interpreted. 

Two related predictions are that for a gap with line nodes that C v ~ H^T 
for H ^ and T <C T c in a polycrystalline sample [108,113] and also that 
the heat capacity scales [114] as a function of H?T. Both of these predictions 
have been confirmed for YBa2Cu307_a; [115] and I^-zSr^CuCU [116]. 

None of the effects described in this section appear to have been studied 
systematically in any of the organic charge transfer salts. Clearly such exper- 
iments could be extremely powerful tools for elucidating the structure of the 
gap in the organic charge transfer salts. 

Thermal conductivity 

In the cuprates, UPt3 and Sr2Ru04 the thermal conductivity is dominated 
by quasiparticles (rather than phonons) over a significant temperature range 
below T c . When quasiparticles dominate the measurements the dependence of 
thermal conductivity of a type II superconductor in the Abrikosov phase on 
the orientation of the magnetic field can yield information about the struc- 
ture of the superconducting order parameter. For a superconductor with an 
isotropic gap one expects that the thermal conductivity, k varies as cos 9 
where is the angle between the heat current and the magnetic field [117]. 
This prediction has been confirmed in Nb [118]. With a thermal gradient along 
the a- [53] or 6-axes [119] of YBa2Cu307_ 2: an additional four fold variation 
in k(9) is observed when a magnetic field is rotated in the basal plane. The 
thermal conductivity is maximal at 9 — 7r/4 + nir /2 for ngZ and therefore 
consistent with l d x 2_ y 2^ superconductivity, i.e., the thermal conductivity is 
maximal with the field aligned in the antinodal direction [120]. In contrast in 
the (low temperature, low field) B phase of UPt3 only the two fold variation 
of k{9) expected for an isotropic gap is observed when the field is rotated in 
the basal plane [121]. This is consistent with either an Ez u or an E± g hybrid 
gap structure [106]. 

In Sr2Ru04 a four fold anisotropy is observed in k(6) as the field is rotated 
in the basal plane [122]. However, this is believed to be a result of the tetrag- 
onal crystal structure rather than a reflection of nodes (or anisotropies) in the 
gap [107]. In particular in Sr 2 Ru04 the fourfold anisotropy is a much weaker 
than the effect in YBa 2 Cu307_ a; , and 20 times small than that predicted [123] 
for a gap with vertical line nodes. 

Izawa et al. have measured the thermal conductivity of k-(ET)2Cu(NCS)2 
as a magnetic field is rotated in the conducting (b-c) plane [93]. They observe a 
small (0.2% for T < 0.6 K) fourfold variation in n{9). This fourfold anisotropy 
has its maximum when the field is at 45° to the b and c axes (i.e., along the 
x and y axes as defined in figure 1.1). Therefore, Izawa et al. propose that 
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the superconducting order parameter of k-(ET) 2 Cu(NCS)2 has t d 3i y symme- 
try. However, one should be cautious of this interpretation because a l d X y 
order parameter transforms like the A g representation of C-zh (see table 1.4). 
Therefore, the anisotropy Izawa et al. observe is has the same symmetry as 
the lattice. Therefore, the observed anisotropy is compatible with variations 
caused by the crystal lattice itself and the observed effect is rather small. In- 
deed Izawa et al. also observed a large twofold anisotropy, which they claim 
is "mainly due to phonons" . 

Furthermore, Izawa et aVs interpretation is based on the assumption that 
there is strong interlayer coupling, which may not correct for k-(ET)2Cu- 
(NCS)a. In particular the theory used to interpret these experiments requires 
that the superconductivity may be described by a (highly anisotropic) three 
dimensional Ginzburg-Landau theory. However, in k-(ET)2Cu(NCS)2 there is 
a very real possibility that the layers are Josephson coupled; in which case 
a three dimensional Ginzburg-Landau theory could not be applied to these 
results. In particular Izawa et aVs interpretation of the data requires that 
vortices are formed when the field is parallel to the highly conducting layer. 

The thermal conductivity of A-(BETS) 2 GaCl 4 and K-(BETS) 2 FeBr 4 has 
also been measured by Tanatar and co-workers [124]. However, in this work 
the field was not rotated in the conducting plane, so we will not discuss it 
further. 

Ultrasonic attenuation 

Because the velocity of sound is much less than the velocity of an electron at 
the Fermi surface (v s <C vp), longitudinally polarised ultrasound is only atten- 
uated effectively by electrons moving almost perpendicular to the direction of 
sound propagation [125]. This makes longitudinal ultrasound a powerful probe 
of the anisotropy of the superconducting gap, indeed longitudinal ultrasonic 
attenuation experiments were among the first to probe gap anisotropy in con- 
ventional superconductors such as Sn [80] . Transverse ultrasound experiments 
arc extremely sensitive to gap anisotropics as the attenuation depends on both 
the direction of sound propagation and the direction of the polarisation. 11 
Such experiments have been important for determining the gap structure of 
both UPt 3 [126,127] and Sr 2 Ru0 4 [128]. 

Measurements of attenuation of longitudinal ultrasound perpendicular to 
the highly conducting planes have been made in k-(ET)2Cu(NCS)2 [59,129, 
130], K-(ET) 2 Cu[N(CN) 2 ]Br [59,131] and k-(ET) 2 Cu[N(CN) 2 ]C1 [60]. A clear 
coupling to the electronic degrees of freedom is seen as the crossover from the 
'bad-metal' to the Fermi liquid is observed at temperatures corresponding to 
the maximum in the resistivity. Indeed the electronic pressure-temperature 



11 For a clear explanation of this phenomena see [107] and [126] and references 
therein. 
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phase diagram of k-(ET) 2 Cu[N(CN) 2 ]C1 has been mapped out using ultra- 
sound [60]. However, none of these experiments were designed to measure the 
anisotropy of the superconducting order parameter. 

Ultrasound experiments arc complicated by the small size and the shape of 
single crystals of organic charge transfer salts [59, 130]. However, for example, 
the crystal of k-(ET)2Cu(NCS)2 studied by Simizu et al. [130] (which was 
reported to be 1.35 x 5.06 x 1.62 mm 3 ) is rather similar size an shape to the 
crystal of UPt3 studied by Ellman et al. [126] (1 x 1 x 2.7 mm 3 ). Therefore, 
there does not seem to be any intrinsic reason why transverse ultrasound could 
not be used to study the order parameter of organic superconductors. Given 
the importance of such experiments [80,126-128] in UPt3, Sr 2 Ru04 and Sn 
transverse ultrasound appear to be an excellent tool to probe the structure of 
the gap in organic charge transfer salts. 

1.5.3 Superfluid stiffness (evidence for a non-BCS groundstate and 
pairing mechanism) 

The defining property of a superconductor is the Meissner effect, i.e., perfect 
diamagnctism [104]. The strength of the Meissner effect is measured by the 
superfluid stiffness, D s = c 2 /47rA 2 , where A is the penetration depth, as the 
superfluid stiffness is the constant of proportionality between the vector po- 
tential of an applied magnetic field and the induced supercurrent. Thus the 
smaller the superfluid stiffness the 'weaker' the superconductivity. In the un- 
derdoped cuprates it is found that T c cx D s : this is the Uemura relation [132]. 

A number of exotic mechanisms have been proposed to explain the Uemura 
relation and it is natural in both the preformed pairs scenario [133] and the 
resonating valence bond (RVB) theory [134]. However, the Uemura relation in 
the undcrdoped cuprates is even predicted within BCS theory [135]. In BCS 
theory D s oc n s /m*, where n s is the density of electrons in the superfluid 
condensate. However, we stress that there is no way to directly measure n s 
and, in general, n s is not required to be the same as the electron density, n, 
even at T = (for example only about 10% of the atoms form the condensate 
in 4 He at absolute zero [136, 137]). However, BCS theory is a weak coupling 
theory and n s = n at T = 0. Recall that the doping, x, is an implicit parameter 
when one compares T c and D s in the undcrdoped cuprates. Therefore, the 
Uemura relation can be derived from BCS theory by noting that T c ~ x, 
which is observed experimentally for small x, and that D s cx n s = n oc x 
which assumes m* is independent of, or only weakly dependent on, doping. 

Pratt and coworkers [138] have systematically measured the critical tem- 
peratures and penetration depths of a large number of organic charge transfer 
salts. This data, along with that of several other groups [139, 140] is shown in 
figure 1.6. It is important realise that the implicit parameters in this plot are 
chemical substitution and pressure, which do not change the filling factor: all 
of these materials are half-filled [141]. The role of pressure, be it 'chemical' 




Fig. 1.6. The small superfluid stiffness observed in organic charge transfer salts with 
low critical temperatures provides clear evidence that the superconductivity is not 
described by BCS theory. The predictions of the phase fluctuation theory proposed 
for the cuprates by Emery and Kivelson [133] is also to overestimate the critical 
temperature. This data is extremely surprising as the data on the left corresponds 
to materials where the interactions are stronger, e.g., the effect mass is larger [141]. 
Thus the superfluid stiffness D s cx 1/A 2 , where A is the zero temperature penetration 
depth, gets smaller (and hence less like the BCS prediction) as we move away from 
the Mott transition. This is unexpected as the spectral weight in the quasiparticle 
peak increases as the correlations become weaker, so we might expect the superfluid 
stiffness to be smallest near the Mott transition [30, 142]. Note that all of these 
materials are half-filled [141] so doping effects cannot be invoked to explain the 
data. This figure is modified from [141], while the data is from Pratt et al. [138], 
Lang et al. [139] and Larkin et al. [140] (colour of citation matches colour of data 
points). The figure shows data for K-ET 2 Cu[N(CN) 2 ]Br (+), k>ET 2 Cu(NCS) 2 at 
several pressures (•), A-BETS 2 GaCl 4 (■), /3-ET 2 IBr 2 (♦) a-ET 2 NH 4 Hg(NCS) 4 (A) 
and K-BETS 2 GaCl 4 (★). 



or hydrostatic, in these materials is to drive them away from the Mott tran- 
sition [28]. Gcnerically, it is also found that increasing the pressure lowers T c . 
Therefore, the data points on the left hand side of figure 1.6 are closer to the 
Mott transition than those on the right hand side. 

As there is no doping of the system BCS theory (and other weak coupling 
theories) does not predict any change in n s as pressure varies T c (whence the 
vertical line in figure 1.6). However, one might object that as we lower pressure 
and drive towards the Mott transition m* increases. If one includes this effect 
the prediction is that the superfluid stiffness is smallest closest to the Mott 
transition. This is the opposite behaviour to that seen experimentally. It seems 
extremely unlikely that using Eliashberg theory to account for strong coupling 
phononic effects will rectify this essential disagreement with experiment [141]. 
Therefore, Pratt et al.'s data [138] provides the clearest evidence (i) for a 
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non-phononic pairing mechanism and (ii) that weakly correlated theories are 
insufficient to explain the observed superconductivity in the organic charge 
transfer salts [141]. 

1.6 Strongly correlated models of superconductivity: 
frustration and RVB 

Most of the early work on superconductivity in the organic charge trans- 
fer salts took weakly correlated approaches. A comprehensive review of this 
work was recently published by Kuroki [143], and, rather than duplicating 
that effort here, we limit ourselves to a few general comments. The two most 
studied weakly correlated theories of superconductivity in the organic charge 
transfer salts are the Eliashberg phononic pairing mechanism and the spin- 
fluctuation pairing mechanism [with calculations most often performed within 
the fluctuation-exchange (FLEX) approximation] . Neither of these approaches 
capture the Mott transition or the large effective mass enhancement seen in 
the organic charge transfer salts. Therefore, one must clearly go beyond a 
weakly correlated description of the superconductivity in order to provide a 
complete description of the full range of behaviours observed in the organic 
charge transfer salts (c.f. sections 1.3 and 1.4). Further, the small supcrfluid 
stiffness observed in the low T c materials (figure 1.6 and section 1.5.3) cannot 
be accounted for within weakly correlated theories [141]. 

Recently a number of groups have proposed strongly correlated theories 
of the organic charge transfer salts. These are based on the Hubbard model 
described in section 1.2.2. A number different methods have been discussed 
[30, 144-149] with rather similar results. Here we focus on the simplest of 
these theories, the resonating valence bond (RVB) or gossamer superconductor 
theory (the two names have been used interchangeably in the literature). 

The celebrated BCS wavefunction [4] is 



where |0) is the vacuum state and the Uk and the t>k are variational parameters. 
The RVB wavefunction is a projected BCS wavefunction, 




(1.11) 



k 



\RVB) = Pq\BCS), 



(1.12) 



where 




(1.13) 



is the partial Gutzwillcr projector. For a = 1 the Gutzwillcr projector removes 
all double occupation from the wavefunction. The a = 1 RVB state has been 
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studied since the early days of high-T c [10, 134]. However, as no double occu- 
pation can occur for a = 1 this state is always insulating for half-filled systems 
and therefore will not give the correct description of the organic charge trans- 
fer salts. However, if we treat a as a variational parameter [150, 151] then 
some double occupation is allowed, and superconducting and metallic states 
may occur. 

The simplest approach to this theory is to make the Gutzwillcr approx- 
imation (as well as the Hartree-Fock-Gorkov approximation implicit in the 
BCS state) in which one enforces the constraints on the fraction of doubly 
occupied sites only on average. This allows one to derive a strongly correlated 
mean-field theory with only a few parameters [150, 151]. Studies of this the- 
ory [30, 144, 146] have shown that the superconducting state: (i) undergoes a 
first-order Mott transition as is seen experimentally (in the mean-field theory 
this is of the Brinkman-Ricc type [152]); (ii) has quasiparticles with a large ef- 
fective masses which are strongly enhanced near the Mott transition as is seen 
in experiments; and (iii) has a strongly suppressed supcrfluid stiffness. The 
RVB theory also predicts that there is a pscudogap above the superconducting 
state in agreement with NMR experiments (see section 1.8.2). Furthermore, it 
has already been shown that the insulating state of the RVB theory supports 
both Neel ordered states [146, 148, 149] like that observed in the insulating 
state of k-(ET) 2 Cu[N(CN) 2 ]C1 and spin-liquid states [30,145,149] like those 
observed in the insulating state of k-(ET) 2 Cu2(CN) 3 and /3'-[Pd(dmit) 2 ]2^ 
(see section 1.8.5). However, to date there as been relatively little work [149] 
on the comparative stability of these states as the frustration is varied. An- 
other important task for these theories yet to be reported is to give a detailed 
explanation of the small supcrfluid stiffness seen in the low-T c materials (see 
figure 1.6 and section 1.5.3). 

1.7 The phase diagram of the Hubbard model on the 
anisotropic triangular lattice 

Eight years ago [25] one of us proposed a speculative zero temperature phase 
diagram for the Hubbard model on an anisotropic triangular lattice. We now 
know significantly more about the phase diagram because there have been 
many studies of this model utilising many different approximation schemes 
and numerical techniques. In figure 1.7 wc sketch a schematic phase diagram 
which summarises current knowledge. 

1.7.1 The Mott transition 

Only a few points on the phase diagram are known exactly. Because of perfect 
nesting the Mott metal-insulator transition occurs at an infinitesimal U for 
both the square lattice (t' = 0) and isolated ID chains (t = 0) [153]. For the 
isotropic triangular (or more correctly, hexagonal) lattice (f = t) there is no 
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nesting at half filling and so the Mott transition occurs at a finite U. The 
critical U has been estimated by a variety of methods [11,30,154,155] and 
seems to be about U = lOt — I5t. We therefore sketch the Mott transition as 
passing smoothly between these three points in our phase diagram. Several 
approximations suggest that the Mott transition is first order throughout the 
(zero temperature) phase diagram, however the recent work of Imada et al. 
[71,72] and the critical exponents measured by Kawagwa et al. [70] has called 
this into question. In particular Imada et al. propose that for some frustrations 
the transition is first order while for others there is a (marginal) quantum 
critical point. All of the other phase transitions shown in the schematic phase 
diagram (figure 1.7) are thought to the second order. 

1.7.2 The superconducting states 

A number of studies based on strongly correlated theories [30, 144-146, 148, 
149] suggest that superconductivity is realised on the metallic side of the Mott 
transition. This superconductivity is mediated by antiferromagnetic interac- 
tions which occur because of supcrcxchangc. As the superexchange interac- 
tions arc only relevant in the large U limit, for small U this interaction will 
vanish and so will the superconductivity, leaving a metallic state for W ^> U. 

The symmetry of the anisotropic triangular lattice is represented by the 
Civ point group, which is summarised in table 1.5. Various calculation meth- 
ods indicate that in the small t' ft limit the superconducting order parameter 
transforms like the A 2 representation of C<i v [24, 144]. This is colloquially re- 
ferred to as l d x 2_ y 2 ' superconductivity (see section 1.5.1 for a discussion of this 
nomenclature). FLEX [143] and weak-coupling renormalisation-group [167] 
studies also predict that the superconducting state has this symmetry. How- 
ever, for t = t' the model has Ce v symmetry. Here 'd x 2_ y 2' state belongs to 
the two-dimensional E2 representation. This has recently led to the proposal 
that for t ~ t' a 'd + id' state (which is a mixed state transforming like a 
complex combination of both the A\ and Ai representations of C 2v ) is re- 
alised [24, 144]. These arguments also suggest that in the large t'/t limit the 
state transforming like the A\ representation of C^v, which may be colloqui- 
ally referred to as a l d xy + s' state, occurs [24, 144]. We therefore sketch these 
different superconducting states in figure 1.7. 

1.7.3 The insulating states 

We now turn our attention to the insulating state. We begin by considering 
the large U/W limit, in which our Hubbard model reduceds to the Heisenbcrg 
model on the anisotropic triangular lattice. This model contains two param- 
eters, J = At 2 /U and J' = At' 2 /U which occur because of of superexchange: 
a second order perturbation effect. All higher order process can be neglected 
because we are in the large U limit. For small J' / J it is well established that 
the Neel state is realised. Series expansions [156] show that the Neel state is 
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J7(J'+J) 




t'/(t'+t) 

Fig. 1.7. Schematic diagram of the proposed zero temperature phase diagram of 
the Hubbard model on the anisotropic triangular lattice at half filling. We argue 
throughout this review that this model provides a qualitative description of the /3, 
P' , k and A phase organic charge transfer salts. The central feature of the phase 
diagram is the Mott transition. We have argued that superconductivity occurs on 
the metallic side of the Mott transition. It has been proposed that the symmetry 
of the superconducting state varies as the frustration is varied and that the super- 
conducting state breaks time reversal symmetry when the lattice is approximately 
hexagonal (t ~ t') [24, 144]. We label the superconducting states by their symme- 
tries (the colloquial names are given a more formal basis in table 1.5). It should 
be stressed that the proposed superconducting states are strongly correlated and 
RVB like rather than weakly correlated BCS states [30, 144-146, 148, 149]. As the 
correlations, which mediate the superconductivity, are reduced we recover a metallic 
groundstate. For U 2> t, t' the Hubbard model at half filling is insulating and the 
spin degrees of freedom can be described by the Heisenberg model. On the basis 
of the calculations reported in [156-160] we have also argued that the insulating 
states should change from Neel order to spiral ordered states to a spin-liquid as i'/i 
is increased, i.e., as the frustration is varied. Near the Mott transition corrections 
to the Heisenberg model, such as ring exchange are important. This may stabilise 
the spin-liquid [161]. This phase diagram includes all of the phases observed in the 
layered organic charge transfer salts at low temperatures. We have shown elsewhere 
in this review (particularly sections 1.3, 1.4 and 1.6) that the finite temperature 
behaviour of the half-filled layered organic charge transfer salts is also described by 
the Hubbard model on an anisotropic triangular lattice. Hence we conclude that 
that the Hubbard model on the anisotropic triangular lattice provides an excellent 
qualitative description of these materials. 



1 Strong electronic correlations in organic charge transfer salts 



35 



J7(J'+J) 




t7(t'+t) 

Fig. 1.8. The same phase diagram as is shown in figure 1.7 with the labels removed 
(for clarity) and the names of some materials added. The position of name indicates 
our guestimate of where that material sits on the proposed phase diagram, the white 
arrows indicate the proposed effect of pressure. This is very rough and accurately 
determining the details of the parameter values corresponding to specific materials 
at particular pressures is a major challenge (see in particular section 1.8.1). The 
following shorthand is used in the figure: k-C1 is k-(ET) 2 Cu[N(CN) 2 ]C1, k-CN 3 is 
k>(ET) 2 Cu 2 (CN) 3 , K-Br is K -(ET) 2 Cu[N(CN) 2 ]Br, k-NCS is ac-(ET) 2 Cu(NCS) 2 , /3- 
I3 is /3-(ET) 2 i3 and /?'-dmit is /3'-[Pd(dmit) 2 ] 2 Z. La 2 Cu0 2 , the parent compound of 
the cuprates, is shown with a large U and on the square lattice (t' = 0) to emphasise 
the similarities between the organics and the cuprates [162] . It is known experimen- 
tally that J' /J m 3 in Cs 2 CuCLt [163]. Cs 2 CuCi4 has a spirally ordered ground 
state at very low temperatures, but shows a region of spin-liquid like behaviour at 
higher temperatures [164] which may indicate that the spin-liquid state is energet- 
ically close to the ground state and that coupling in the third dimension plays an 
important role in the real material. k>(ET) 2 Cu 2 (CN)3 shows a spin-liquid ground 
state [165], while the /3'-[Pd(dmit) 2 ] 2 Z salts seem to be on the boarder between 
magnetic ordering and spin-liquid behaviour. The nature of this magnetic ordering 
is not yet clear. Pressure drives k-(ET) 2 Cu[N(CN) 2 ]C1 from a the Neel state via 
the superconducting state to normal metal. While K-(ET) 2 Cu[N(CN) 2 ]Br is so close 
to the Mott transition that it can be driven insulating by deuteration of the ET 
molecule. Huckel calculations suggest that /3-(ET) 2 i3 (and other f3 phase salts) have 
t' > t [166]. If our speculative phase diagram is correct then this suggests that they 
have a different pairing symmetry to the K-phase materials [24, 144] . This also sug- 
gests that the superconducting states of k-(ET) 2 Cu 2 (CN)3 and /3'-[Pd(dmit) 2 ] 2 Z 
may break time reversal symmetry [24, 144] . 
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line 
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Table 1.5. The character table, symmetry required nodes and some basis functions 
of the irreducible representations of the point group C 2v . C 2v represents the sym- 
metry of the anisotropic triangular lattice. The functions lk, X k and Yk may be 
any functions which transform, respectively, as 1, t and k y under the operations of 
the group and satisfy translational symmetry. The operations in the group are the 
identity, E, rotation about the z-axis by n, C 2 , and reflection through the planes 
kx = ik y , a ± . A brief explanation of characters is given in the caption to table 1.3. 
Note that inversion symmetry is not an operation of the group. This might suggest 
that singlet and triplet states are not differentiated by symmetry. However, as we 
are discussing the symmetry of a two dimensional lattice, rotation by n (which is 
an operation in the group) differentiates between singlet and triplet states. 



stable for J' / J < 0.7. For J = we have uncoupled one-dimensional chains 
along the diagonal on the unit cell. The Heisenberg model can be solved ex- 
actly in one-dimension and the ground state is a spin-liquid with deconfincd 
spinon excitations [168]. The question of what happens when there is a finite 
coupling to a second dimension (finite J) is extremely subtle. However, recent 
work [147, 156, 169] suggests that this physics docs survive to a significant 
degree in frustrated systems and that a spin-liquid occurs in some parts of 
the phase diagram with J' > J. Therefore, we label the large J'/ J region of 
the phase diagram 'qlD?'. 

Various theoretical studies of the Heisenberg model on the anisotropic 
triangular lattice have been performed. The methods used include linear spin 
wave theory [157] and large- AT mean field theory [158]. Both of these methods 
indicate that when frustration destroys Neel order [which has the ordering 
wavevector q = (tt, tt)] a spiral state [with ordering wavevector q = (q, q)] 
occurs. In the spiral state q varies from tt, at the critical frustration where the 
Neel state gives way to spiral order, to ir/2 at large frustration which is the 
classical value for uncoupled chains. At J' = J these theories give q = 3n/2 
which describes the 'classical 120°-state': the classical solution of the isotropic 
triangular lattice [159]. In the classical 120°-state the spins on neighbouring 
sites align at an angle 120° apart from their nearest neighbours. It is widely 
believed that the quantum analogue of T20°-state' is the true ground state of 
the Heisenberg model on the isotropic triangular lattice [156, 160]. By analogy 
it may be argued that quantum-spiral states are the ground states of the 
Heisenberg model on the anisotropic triangular lattice in the regime J' ~ J. 
This is certainly what is suggested by the large- studies [158]. However, it 
is not yet clear how or at what parameter values the spiral state changes into 
the spin-liquid. 
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At lower values of U/W (i.e., near the Mott transition) higher order per- 
turbation processes cannot be neglected. For example, ring exchange processes 
(which is of order t 4 /U 3 ) [161] and hopping around triangular clusters (which 
is of order t 2 t' /U 2 ) [11] may become important. Both of these processes favour 
RVB spin-liquid states. This may explain the observation of a spin-liquid state 
in insulating phase of k-(ET) 2 Cu 2 (CN) 3 [165] and /3'-[Pd(dmit) 2 ] 2 Z [170] (see 
section 1.8.5). Therefore, we add an RVB spin- liquid phase to our sketch of 
the phase diagram intervening between the quasi-one-dimensional phase and 
the spiral order. We also include the enhancement of the region of stability of 
this phase near to the Mott transition. 

1.8 Some outstanding problems 

Before drawing our conclusions, we briefly indulge ourselves by discussing 
what we see as some of the major problems in the field that we have not 
discussed significantly above. 

1.8.1 Deducing parameters for the minimal models from first 
principles electronic structure calculations 

We have argued above that the physics of the organic charge transfer salts can 
be understood in terms of the Hubbard model on the anisotropic triangular 
lattice. Clearly it is vitally important, for this program, to know accurately 
what values of the parameters of this model (t, t' and U) correspond to a 
given material at a given pressure. Understandably, given the chemical com- 
plexity and large unit cells of the organic charge transfer salts, most theoretical 
work on the band structure of these materials has used the Huckel model [15]. 
(There has also been considerable effort expended to determine the values of t 
and t' experimentally [171].) However, recent advances in computational speed 
and the computational efficiency of electronic structure codes [172] have al- 
lowed the first density functional theory (DFT) studies of the band structures 
of organic charge transfer salts to be performed [17-20]. The calculations are 
significantly more accurate than the Huckel calculations and there is a great 
need for systematic studies of the band structures of a range organic charge 
transfer salts. In particular a detailed mapping from the experimental param- 
eters of 'chemical' and hydrostatic pressure to the theoretical parameters of t 
and t' would enormously benefit the field. Further, this would greatly improve 
our ability to perform quantitative tests of theories of the Hubbard model in 
the organic charge transfer salts. 

To complete the mapping between the experimental and theoretical pa- 
rameters we also need to know how the dimer U varies with 'chemical' and 
hydrostatic pressure. This is a much more difficult task. To date the most com- 
mon approach has been to estimate U from the intra-dimer hopping integral. 
We have shown in section 1.2.2, that this is not only inaccurate (as it is based 
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on extend Huckel calculations) but also incorrect as correlation effects, not 
the intra-dimer hopping integrals, determine the dimer U. Therefore, more 
accurate calculations are required. The 'bare' Hubbard-C/ of a dimer in vacuo 
can be straightforwardly determined form DFT as it is the second derivative 
of energy with respect to charge (U = d 2 E/dq 2 ). However, in a crystal the 
effective Hubbard-[7 is greatly decreased as the second derivative of energy 
with respect to charge contains a large contribution due to the polarisability 
of the crystal. The problem of calculating the reduction in the Hubbard-[7 
due to the polarisability of the molecules is simply a self-consistent problem 
in classical electrodynamics. This approach has been successfully applied to 
alkali doped fullerides [47, 173]. Thus the large size and low electron density 
of molecular crystals is a singular advantage when it comes to calculating the 
Hubbard-?/. Indeed this is such an advantage that, we believe, this makes 
molecular crystals an ideal class of materials in which to study first principles 
approaches to strongly correlated effects. We believe that the time is now ripe 
to apply this method to the organic charge transfer salts and thus to map 
from 'chemical' and hydrostatic pressure to the Hubbard-t/. 

If both of these procedures were accurately carried out then we would have 
a direct mapping between the Hubbard model and experiment. This would 
allow very direct testing of the hypothesis that the Hubbard model provides 
the correct microscopic description of the organic charge transfer salts. It 
would also allow for the experimental testing of the various approximations 
used to study the Hubbard model. Therefore, this mapping would not only 
be extremely important to the community interested in the organic charge 
transfer salts for their own sake, but would make the organic charge trans- 
fer salts an even more important test-bed for theories of strongly correlated 
quantum many-body systems than they are already. In figure 1.8 we sketch 
our best guess of which parameters correspond to which materials, based on 
the calculations performed thus far and comparison of the observed behaviour 
of these materials with figure 1.7. Given the comments above this is clearly a 
very rough procedure. 

1.8.2 Is there a pseudogap? 

NMR experiments on the organic charge transfer salts show a large decrease in 
the spin lattice relaxation rate, 1/T±T, the Knight shift, K s , and the Korringa 
ratio, IjTxTKl below about 50 K [83-85] (see figure 1.9). In this section we 
ask whether these experiments indicate the opening of a gap-like-structure at 
the Fermi energy. The pseudogap phase in the cuprates has attracted a great 
deal of attention [174], yet there has been very little work on understanding 
the origin of these experimental effects in the the organic charge transfer salts. 
Clearly a good starting point to investigate this effect in the organic charge 
transfer salts would be to follow a number of the experimental and theoretical 
approaches that have proved fruitful in investigating the pseudogap in the 
cuprates. 
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Fig. 1.9. Is the suppression of the nuclear spin relaxation rate, 1/Ti, evidence for 
a pseudogap? This data (from [175]) shows the temperature dependence of 1/T\T 
in K-(ET) 2 Cu[N(CN) 2 ]Br at various pressures [1 bar (•), 1.5 kbar (A), 3 kbar (o) 
and 4 kbar (o)]. Rather similar effects are seen across a range layered organic charge 
transfer salts [85]. In materials that are close to the Mott transition a peak in 
1/T\T is observed in at around 50 K. This is not expected in weakly interacting 
systems. However, a pseudogap is predicted by the RVB theory of superconductivity 
[30] and appears to be a rather natural feature of extensions to dynamical mean 
field theory which include some spatial correlations [176,177]. A detailed theoretical 
understanding of this data is still lacking, as are further experiments probing the 
possibility that there is a pseudogap in the organic charge transfer salts. 



One of the most important questions about the pseudogap in the cuprates 
has been whether, and if so how, it relates to superconductivity. The Nernst 
effect is much larger in a type-II superconductor than it is in a normal metal 
(because vortices carry entropy extremely efficiently). Therefore, the discov- 
ery of a large Nernst effect in the pseudogap regime [178] of the cuprates is 
suggestive of superconducting fluctuations playing an important role in the 
pseudogap [133]. In this context studies of the temperature dependence of the 
Nernst effect in materials which show the unexpected NMR behaviour {e.g., 
/c-(ET) 2 Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 } would be extremely interest- 
ing. 

Theoretical avenues could also be productively pursued. In particular we 
need to discover the correct phenomenological description of the NMR data. 
The first question which needs to be addressed is can the NMR data be ex- 
plained without invoking a gap in the density of states? For example, the 
maximum in \jT\T occurs at a temperature very close to that where the 
crossover from the 'bad- metal' to the Fermi liquid occurs. Could the origin of 
the maximum then be simply related to recovery of a Fermi-liquid Korringa- 
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like behaviour from the local moment (Hcisenbcrg) physics associated with 
the 'bad-metal' [179-181] (charge fluctuations happen on significantly slower 
time-scales than spin fluctuations in the 'bad-metal')? Or, is a gap in the den- 
sity of states necessary to describe the observed phenomena? A phenomeno- 
logical description of the data will obviously provide clues to the appropriate 
microscopic description. Given the success of DMFT in describing many of 
properties of the organic charge transfer salts (see section 1.3.1), extensions to 
DMFT to include non-local effects seem an obvious avenue to explore [176]. 
Cellular DMFT (CDMFT) [176] and the dynamical cluster approximation 
(DCA) [176] are two such approaches. These approaches have already pro- 
vided significant insights into the pseudogap in the cuprates, in particular 
the so-called Fermi arcs seen in angle resolved photoemission spectroscopy 
(ARPES) experiments [182,183] on the cuprates appear to be a natural feature 
in CDMFT [177]. Therefore, studies of the Hubbard model on an anisotropic 
triangular lattice within CDMFT or the DCA would be an extremely inter- 
esting approach. 12 It is interesting to note that many of the features present 
in the CDMFT description of the pseudogap are captured (in a somewhat 
cruder form) by the RVB theory (discussed in section 1.6). In particular the 
RVB theory predicts a pseudogap [30,144,146,151]. Therefore, a detailed un- 
derstanding of the NMR experiments on the organic charge transfer salts will 
provide a stringent test of this theory. 

1.8.3 Low T c materials 

Most of the attention to superconducting has been focused on materials with 
high transition temperatures and the region around the Mott transition. But, 
recently, the work Pratt, Blundcll and coworkers [138,141,142] has shown the 
the importance of materials with low critical temperatures for understanding 
superconductivity in the organic charge transfer salts (see figure 1.6 and sec- 
tion 1.5.3 for a discussion of this work). However, there remains very little data 
on the low-T c materials. We believe that an extensive understanding of the 
phenomenology of these low-T c materials is a vital prerequisite of a more de- 
tailed understanding of the superconducting states of all of the organic charge 
transfer salts. Therefore, there is a desperate need for more experimental and 
theoretical work on these materials. 

1.8.4 Synthetic chemistry 

It goes, almost, without saying that condensed matter physics cannot proceed 
without high quality samples. As such, synthetic chemistry plays the central 
role in the study of organic charge transfer salts, for, without it, the field could 
neither exist nor progress. Historically, much of the synthetic work in this 

12 One CDMFT study of the Hubbard model on an anisotropic triangular lattice 
has appeared in the literature but this problem was not addressed in that paper [149] . 
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field has focused on discovering new organic charge transfer salts. However, 
we would like to stress the importance of the less obviously glamorous work 
of producing higher quality samples of materials which have already been 
extensively studied. As we have seen in this review there are many outstanding 
problems in the flavours of organic charge transfer salts that have been with us 
for the last twenty years. Many of these problems would be made significantly 
more tractable if new kinds of data were available. For example, a tricrystallinc 
experiment, of the type so important in the cuprates [98], could lay to rest 
once and for all questions about the pairing symmetry in the organic charge 
transfer salts. 

Perhaps the most pressing need for improvement in sample quality is size. 
Neutron scattering has been extremely important for our understanding of 
the cuprates [184], Na^CoC^ [11, 185] and many other strongly correlated 
materials. This is mostly because neutrons provide direct information about 
the spin degrees of freedom which are of vital importance in strongly correlated 
materials such as the organic charge transfer salts. Therefore, perhaps the 
single greatest impediment to the field is the lack of high quality single crystals 
for neutron studies. Taniguchi and coworkers have recently demonstrated that 
large single crystals can be grown [186] and one can only eagerly await the 
results of neutron scattering experiments on this sample and the growth of 
similarly large samples of other salts. 

1.8.5 Do k-(ET) 2 Cu 2 (CN) 3 and/or ^'-[Pd(dmit) 2 ] 2 Z have 
spin-liquid ground states? 

A recent series of experiments on k-(ET)2Cu2(CN)3 [165, 187] by Shimizu 
and coworkers has sparked considerable interest in the theoretical community 
[30,56,144,145,149,161,188-195]. This is because Shimizu et al. have gathered 
strong evidence that k-(ET)2Cu2(CN)3 may have a spin-liquid ground state. 
A spin-liquid is a state in which there exist well formed local moments which 
do not magnetically order. Importantly a spin-liquid should posses all of the 
symmetries of the crystal. In a particularly beautiful demonstration of the 
spin-liquid behaviour in k-(ET) 2 Cu2(CN)3 Shimizu et al. compared the spin 
susceptibility of k-(ET) 2 Cu 2 (CN) 3 with that of k-(ET) 2 Cu[N(CN) 2 ]C1 (see 
figure 1.10). Both k-(ET) 2 Cu 2 (CN) 3 and k-(ET) 2 Cu[N(CN) 2 ]C1 have well 
formed local moments with spin-i. In both systems the Heisenberg exchange 
interaction is about 250 K [25,165]. However, the frustration is significantly 
larger in /c-(ET) 2 Cu 2 (CN) 3 than it is in k-(ET) 2 Cu[N(CN) 2 ]C1 [56,165,196]. 
The susceptibility of k-(ET) 2 Cu[N(CN) 2 ]C1 diverges at -25 K as the material 
undergoes the Neel transition. In contrast no magnetic transition is observed 
in k-(ET)2Cu2(CN)3 down to 32 mK (the lowest temperature studied by 
Shimizu et al.). 

Rather similar behaviour (also shown in figure 1.10) has been observed 
by Kato et al. [170] in a the series /3'-[Pd(dmit) 2 ]22', where Pd(dmit)2 is 
the acceptor molecule l,3-dithiol-2-thione-4,5-dithiolate (C3S5), Z=MjzaY or 
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Fig. 1.10. Evidence for spin-liquid ground states in k-(ET)2Cu2(CN)3 and ff- 
[Pd(dmit)2]2-^- The left panel (modified from [165]) compares the temperature de- 
pendencies of the magnetic susceptibilities, %, of k-(ET)2Cu2(CN)3 and k-(ET)2- 
Cu[N(CN)2]C1. Both measurements were performed in the ambient pressure Mott 
insulating phases of the two compounds. In k-(ET)2Cu[N(CN)2]C1 a rapid upturn in 
the susceptibly is observed at ~25 K as a result of the Neel transition. Remarkably 
no magnetic ordering transition is observed in k-(ET)2Cu2(CN)3 down to the lowest 
temperatures studied (32 mK). Fits of series expansion data (lines in figure) to the 
measured susceptibility [165] show that the exchange interaction, J, is about 250 K 
in both materials. This had led to the suggestion that it is the greater frustration 
in k-(ET) 2 Cu2(CN) 3 that leads to a spin-liquid ground state. [In k-(ET)2Cu 2 (CN) 3 
J' I J ~ 1 whereas in k-(ET) 2 Cu[N(CN) 2 ]C1 J'/J < 1.] The right panel shows sim- 
ilar measurements in /3'-[Pd(dmit)2]2^ [170]. This data is also well described by 
series expansions for the isotropic triangular lattice (solid lines in figure with J 
marked). Except for X=Et2Me2Sb the /3'-[Pd(dmit)2]2^ salts do eventually order 
magnetically, but this suggests that it may be possible to use chemistry and/or uni- 
axial stress tune the /3'-[Pd(dmit)2]2^ salts across the antiferromagnetic/spin- liquid 
phase boundary. Understanding the microscopic nature of these spin-liquid states 
and why they are stabilised over the 120°-state which appears to be the ground state 
of Heisenberg model on the isotropic triangular lattice present major challenges to 
theory. 



Et 2 Mc 2 F, Y=P or Sb, Mc=CH 3 and Et=C 2 H 5 . So far, this has not attracted 
so much attention from theorists as the k-(ET)2Cu2(CN)3 results, but these 
materials are also excellent candidates for the observation of spin-liquid states. 
The susceptibility of the /3'-[Pd(dmit)2]2^ compounds is much more reminis- 
cent of that of k-(ET) 2 Cu 2 (CN) 3 than that of k-(ET) 2 Cu[N(CN) 2 ]C1. In par- 
ticular the Et2Me2Sb salt shows no indications of magnetic ordering down to 
4.3 K (the lowest temperature studied thus far, although magnetic impurities 
cause a Curie tail and make this somewhat difficult to substantiate) [197]. 
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Further, the range of Curie temperatures seen in the salts that do eventually 
order magnetically [197] suggests that varying the cation, Z 1 in the /3'-[Pd- 
(dmit^^-Z series allows one to tune the frustration and thus the proximity to 
the spin-liquid state. 

The clear theoretical challenge is to explain why a spin-liquid ground state 
is stable in these materials. 13 Simply appealing to the geometrical frustration 
inherent in the triangular lattice is insufficient, as it is all but certain that 
the 120°-state is the ground state of Heisenberg model on a triangular lattice 
(see section 1.7). A great deal of theoretical effort has already been expended 
on this problem [30,56,145,149,161,188-195]. The idea that proximity to 
the Mott transition allows perturbation terms not included in the Heisenberg 
model to become large seems particularly promising. However, this question 
is far from settled yet and will doubtless be the basis of much debate in the 
future. 



1.9 Conclusions 

We have seen above that the predicted behaviour of the Hubbard model on 
the anisotropic triangular lattice provides good qualitative agreement with 
experiments on the organic charge transfer salts. This success is not just lim- 
ited to the k phases but the model applies equally well to the (3, (3 1 and A 
phases: in spite of the chemical and structural differences between these ma- 
terials the physics is essentially the same. As the exact solution of this model 
is not yet known various approximation schemes must be used to discover the 
true behaviour of the model. Perhaps the most notable successes of this model 
are the explanation of the metallic state and the Mott transition in terms of 
dynamical mean field theory (DMFT) discussed in section 1.3.1. However, 
recent theories of superconductivity based on the resonating valence bond 
(RVB) state appear to explain many of the features of the superconducting 
state. An important test of the RVB theory will be a detailed comparison 
of its prediction of a pseudogap with experiments (such as NMR or Ncrnst 
effect) on the strongly correlated metallic state just above T c . An important 
challenge for this theory is to explain the small superfluid stiffness seen far 
from the Mott transition, i.e., in compounds with low critical temperatures. 

There are many exciting challenges facing the field. Some are old - like 
finding an experiment which decisively settles the questions about the pair- 
ing symmetry; while others are extremely new - like understanding the spin- 
liquid state in «-(ET)2Cu2(CN)3. An important challenge is also to increase 
the quantitative detail of the predictions of theory. This is required, not only 
because of intrinsic interest in the organic charge transfer salts themselves, 
but also because the organic charge transfer salts provide a wonderful test-bed 



Even if the compounds do order at some extremely low temperature the spin- 
liquid state is clearly at least very energetically competitive. 
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for many of the most important ideas in the theory of strongly correlated ma- 
terials. In particular, the Mott transition remains perhaps the phenomena of 
central importance to strongly correlated physics and the organic charge 1 trans- 
fer salts are one of only a handful of systems where the Mott transition can be 
driven by varying U/W as Mott originally envisioned [198] rather than by dop- 
ing. We have highlighted one possible approach to increasingly quantitative 
prediction, that of paramcterising minimal models, although there are several 
other approaches which might be profitable, such as LDA+DMFT [199]. 

We hope that this review has made it clear that the organic charge transfer 
salts are a playground for quantum many-body physics. 14 Many of the most 
important phenomena in strongly correlated physics are found in the organic 
charge transfer salts, for example, the Mott transition, unconventional super- 
conductivity, frustrated antifcrromagnetism and spin-liquids. Further, the or- 
ganic charge transfer salts are exceptionally clean systems (as is evidenced by 
the beautiful quantum oscillations observed at low temperatures [40,41,171]). 
But, perhaps, what makes the organic charge transfer salts most attractive 
is the relative ease with which the strength of correlations can be controlled 
via hydrostatic pressure or chemistry [8,28]. Thus, the study of the organic 
charge transfer salts is an important branch of strongly correlated physics. 
The lessons learned from the organic charge transfer salts are already hav- 
ing significant impact in other strongly correlated systems. For example, the 
spin- liquid phase of CS2CUCI4 has a great deal in common with the insulating 
phases of k-(ET) 2 Cu 2 (CN) 3 and /3'-[Pd(dmit) 2 ] 2 Z [56] and the 'Curie- Weiss 
metal' seen in Na^OoO^ is essentially a doped analogue of the 'bad metal' 
observed in the organic charge transfer salts [11]. All of these phenomena 
arise from the interplay of strong correlations with frustration which can be 
so elegantly studied in the organic charge transfer salts. 
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